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ABSTRACT 


This report is presented in two parts. .The first part describes 
the results of the analytical portion of this program to establish the ex- 
pected degree of thrust misalignment in an electric bombardment ion . . 
thruster. The analysis includes the probable misalignment that might 
occur due to mechanical displacement and time, and its effect on the 
thrust vector. An analysis of mechanical displacement due to thermal 
and mechanical stresses on the screen and accel electrodes is presented, 
as is an analysis of the design of thrust stands. The second part de- 
scribes the hardware portion of this program during which a representa- 
tive thruster and thrust stand were fabricated, assembled, and tested. 

A 30 cm, 2 kW thruster was built with an accel electrode that is capable 
of being mechanically and remotely displaced in the axial, tilted, radial, 
and translational modes. The thrust stand is capable of resolving the 
thrust direction to ±1/4° of angle and torques about the thrust axis of 
• 5 x -1-0 inj. -lb. 
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I. INTRODUCTION AND SUMMARY 


As a resuLt of comprehensive studies over the past four years, 
the fundamental design concepts associated with solar powered electric 
propulsion systems for interplanetary missions have become fairly well 
defined. It is known, for example, that the engine array will involve a 
cluster of thruster modules which may be arranged on a frame which 
can be translated in two directions (or other suitable freedom of motion) 
relative to the spacecraft. This translational motion is necessary in 
order to assure that the center of thrust will pass through the center of 
mass of the spacecraft, for any configuration of active and inactive 
thrusters. In addition to changes in the thrust axis due to switching of 
thrusters, the vector from an individual thruster will not necessarily 
be* along the axis and could change slightly in magnitude and direction 
with time. The effect of these misaligned thrust vector components 
around the pitch and yaw axes can be negated by the translation mech- 
anism mentioned above. The component around the roll axis, however, 
can be compensated only by the spacecraft attitude control system or 
by gimballing some or all of the thrusters. Therefore, it is of consid- 
erable importance to SEP spacecraft design to know the expected thrust 
vector misalignment as a function of various thruster operational param- 
eters, mechanical displacements, and time. 

The objectives of this program were to establish by analysis the 
expected change in thrust vector from an electron bombardment thruster 
as a result of all known or anticipated causes; to build a representative 
thruster which was capable of mechanically producing electrode mis- 
alignments; and to build a thrust stand capable of resolving the thrust 
direction and the torque produced by the misalignment. The analysis 
was devoted to expected electrode displacements, determination of the 
effect of these displacements on the ion trajectories and thereby on the 

/ r 

thrust vector, and the design of apparatus. The building of the hardware 
for the thruster and thrust stand was devoted to the verification that 
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each component part of the thruster, the electrode displacement mech- 
anism, and the thrust stand was capable of performing the tasks for 
which it was intended. 

The misalignment of the thrust vector of an array of thrusters 
can arise from a change in the magnitude of the axial thrust or from the 
transverse component of thrust from a single thruster. In turn, the 
total thrust (axial and transverse) of a thruster is made up of the thrust 
elements resulting from each aperture in the accelerator electrode 
system. 

In the analysis of thrust vector misalignment, a comprehensive 
digital computer program was developed in order to determine the ion 
flow and trajectory conditions for a single hole as a function of the rela- 
tive positions of the accelerator and screen electrode and ion current 
density. In this computer study the upstream plasma- sheath (i. e. , the 
source of ions) received careful attention; the program allowed changes 
in the shape and position of the sheath in the screen aperture. The 
position of the plasma emitting boundary was determined by an iterative 
technique, after which detailed trajectories were run. It was found that 
the sheath position and the average direction of the trajectories varied 
as the electrodes were perturbed in position. For changes in spacing 
between accel and screen electrodes, the principal cause of the change 
in current was the variation in plasma boundary area as this surface 
advances or recedes in the tapered screen aperture. With transverse 
electrode displacement the ion trajectories vary, with relatively little 
distortion of the plasma boundary. These relations, between ion optical 
characteristics and electrode displacement, are expressed by a number 
of linear coefficients. 

Several possible causes of electrode misalignment were con- 
sidered. Of these, the one which is inherent in this type thruster and 
which is of greatest importance is the warping or buckling of the screen 
and accel electrodes due to thermally induced stresses. Other causes 
(e. g. , of transverse or rotational misalignment) will result from 
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assembly errors -or from a shift caused by..vibration; these causes. are 
not amenable to accurate. analysis but must be evaluated accurately by 
experiments. The thermal buckling problem received considerable 
attention. The temperature distribution in the screen and accel elec- 
trodes due to various heat loads, heat input distributions, etc., was 
analyzed. The effect of these temperature distributions on the mechani- 
cal distortion of the electrodes was determined by means of a digital 
computer program. The predicted temperature distribution was com- 
pared with experimental data in one case, with good agreement. There- 
fore, by using the linear coefficients determined from the ion optical - 
study together with expected electrode displacements, the expected 
change in thrust is determined. It is found, for example, that a change 
in axial thrust as high as 18% can be expected as a result of electrode 
warping. 

The program involved the design and fabrication of experimental 
equipment for measuring the thrust vector of an electron bombardment 
thruster. A thruster was designed — of 30 cm, 2 kW size — to be repre- 
sentative of flight-type state-of-the-art mercury bombardment thrusters. 
During the early analytical phase of the program a parallel HAC IR&D 
program allowed the experimental verification of the suitability of a 
similar thruster for the intended purpose. 

During the evaluation of a number of different forms of thrust 
stands capable of measuring the four components of thrust desired, a 
new concept was conceived, analyzed in detail, and adopted for use. 

This thrust stand design is novel in that the platform is floated on a 
pool of mercury with suitable linear force transducers to allow direct 
measurement of all transverse components of thrust. Comprehensive 
analysis showed that it exhibits the required sensitivity and stability. 

The thrust stand was fabricated and the design details of the components 
and system are presented in this report. 
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In this report the Lon optical study is discussed first, followed 
by an analysis of causes of electrode misalignment (principally thermal 
warping). In Section IV these two analyses are combined to yield the 
expected thrust change. Thrust stand design analysis is given in Sec- 
tion V, The thruster and thrust stand design are given in detail in 
Section VI and VII. Finally, thrust stand performance is presented in 
Section VIII. 
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II. 


ION OPTICAL, STUDY 


A . Introduction 

The total thrust from an ion thruster consisting of an array of 
circular apertures is given by the time rate of change of the momentum 
of the ion particles exhausted from it. The velocity, and hence momen- 
tum, of each ion, is determined by the electric fields it experiences 
during its motion. Thus, a perturbation of these fields caused by an 
electrode misalignment will perturb the total thrust (magnitude and 
direction). A change in the total ion current will cause a similar change 
in total thrust. The purpose of the work reported here was to determine 
the exhaust momentum for various electrode perturbations and ion cur- 
rents, and thus allow calculations of the expected variation in the thrust 
ve ctor . 


A number of the basic electrode misalignments considered are 
illustrated in Fig. 1. These were chosen to include the expected types 
of perturbations which an accel electrode can experience in assembly or 
operation of a thruster. 

In Fig. 1(a) the electrode spacing is decreased while the rotational 
symmetry is preserved. This results in an increase in the magnitude of 
the thrust from this aperture, which may cause the direction of the thrust 
vector from the thruster as a whole to rotate. In Fig. 1(b), the electrodes 
are displaced normal to the beam axis, thus destroying its symmetry 
and causing a change in the magnitude and direction of the thrust vector. 
The case illustrated in Fig. 1 (c) can be considered as the sum of the two 
effects above, which alters both the magnitude and direction of thrust 
ve ctor. 


It can be seen that the rotational symmetry of the flow generally 
will be destroyed and the perturbed flow will be formed about a curvilinear 
axis. It is not possible, in general, to determine this type of flow exactly, 
either with the analog or digital computer. However, a first order solu- 
tion for the thrust vector has been formed by a comb inationof the axisymmetric 
and planar solutions obtained in the digital computer. 
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. Fig. i. Basic electrode misalignments. 
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These ion flow calculations involve two fixed (in shape) bound- 
aries (viz., the screen and accel electrodes) and two fluid boundaries 
(viz. , the upstream and downstream plasma surfaces). Fortunately 
it was found that the shape and position of the downstream boundary 
(the edge of the neutralized region) was relatively insensitive to up- 
stream conditions and vice versa; its small changes did not affect the 
ion flow. On the other hand, the upstream or ion source boundary 
exerted a predominant effect on the ion current and trajectories; thus 
considerable effort has been spent in assuring that the technique of 
treatment of this surface is correct. This factor is also important 
in this study in that the consideration of the fluidity of the upstream 
boundary is an important difference between this thrust vector analysis 
and that which has been carried out previously. 

This section will describe first the techniques used in the ion 
optical study; the analysis of the effect on ion flow of the several elec- 
trode perturbations listed above will then be given. This analysis will 
evolve linear coefficients between thrust and electrode displacement, 
and these data will be used in Section IV to calculate the expected thrust 
vector changes. 

B . Basic Techniques Employed 

1 . Digital Computer Program 

A digital computer program which determines the trajec- 
tories of ions, assuming either rotational or planar symmetry, has been 
used in this work. In most problems involving space charge flow, the 
boundary conditions necessary for solution are prescribed on a given 
surface, on which the potential variation or the normal derivative is 
specified. On the other hand, problems involving plasma boundaries 
involve "free 11 surface in the sense that the location of the boundary 
must be found, given the nominal potential variation and the normal 
derivative. This latter class of problem can be solved only by specifying 
a boundary contour and by refining it in an iterative manner until the 
desired boundary conditions are obtained. 



The basic computer program uses finite difference techniques 
to solve Laplace's and Poisson's equations using boundary data speci- 
fied by the user. Two types of boundaries are of interest: a Dirichlet 
boundary where the potential is known, and a Neumann boundary 
where the normal derivative of the potential is known. In order to 
specify these boundary data, the area of the accelerator is overlaid 

with a suitable mesh and the potential at each mesh point is calculated 

2 

in terms of its neighbors. Hence the Laplace equation V V = 0 at 
point R, Z within the problem becomes RV^ i z + Z- 1 

+ < R + */ 2 » V R + 1, Z + < R - l/2 > V R, Z+l ' 4RV R, 2=0 in terms of finite 
differences, where ^ the potential at the mesh point R, Z. The 

ion source is then divided into segments of equal length and the emission 
is determined by Child's law, using the potential as calculated above at 
a given distance from the surface . 

In practice, the region immediately in front of the emitter is known 
to be highly unstable and a ft suppression factor n must be used in the cal- 
culation of the emission to ensure that the solutions converge. The 
trajectory of a characteristic ion from each segment is then determined 
using the electric fields associated with the Laplace solution. Charges 
appropriate to the current flow are then set in at each mesh point. The 
full Poisson solution can then be determined in terms of finite differences. 
The above process is repeated until the solutions converge. 

Since the emitting surface in this case is a plasma, its shape and 
location are determined by the "self- consistent 11 electric fields at its sur- 
face. If the plasma sheath is considered to be thin and the surface well 
defined, its boundary conditions can be defined as 

V (r^) = const = V plasma 

E (r ) = 0 
s 

P (r fi ) = 0 

where V(r g ), E{r g ) and p(r g ) are the potential, electric field, and 

spacecharge, respectively, at r where r refers to points within 

s s 
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the plasma. A further boundary condition to be satisfied at the up- 
stream boundary is that the ion emission from the plasma should be 
equivalent to the self-consistent space-charge flow from its surface. 

In general, iterative methods must be used to determine the 
location and shape of the upstream and downstream boundaries self- 
con si stently with the flow. The method used was to assume a trial 
boundary of the correct potential and then successively deform this 
until the other boundary conditions were satisfied. An electrolytic 
tank in conjunction with an analog computer was used to obtain an ini- 
tial estimate of the upstream boundary shape and position. For this 
work an assumed downstream plasma boundary was used and the method 
of solution was similar to that described in Ref. 1. The boundary thus 
found was then used in the digital computer in order to determine the 
shape and location of the downstream plasma boundary. 

Z . The Unperturbed Accelerator 

In the electron bombardment thruster the ions are ex- 
tracted from a neutral plasma and given a directed velocity by a suit- 
able electrode system. At an appropriate location outside the thruster 
the accelerated ions are mixed with electrons so that a neutral stream 
of charged particles is ejected from the spacecraft. 

The design of the electrode system is governed by several con- 
siderations. The available ions must be focused efficiently so that the 
accelerated stream of ions has little transverse spread. At the same 
time, direct impingement on the accelerating electrode must be kept to 
a minimum to prevent it from being sputtered away. In addition, the 
mixing of electrons must be accomplished so that backstreaming of the 
electrons to the screen electrode is avoided. The accelerating electrode 
system also must be capable of drawing the required current density from 
the ion emitter. 

The transition from the neutral plasma at a uniform potential to 
the region of pure ion flow, where fields exist, occurs within a very 
small dividing region within which the electrons having energy spread 
«kT are reflected back into the plasma. It is therefore convenient to 



regard the plasma as being bounded by a negligibly thin sheath, the 
Debye sheath. Since the plasma is a region of uniform potential, the 
sheath defines an equipotential; by the same definition the electric 
field normal to the sheath must vanish on the side facing the plasma 
and also on the side where the ion flow takes place. 

The sheath thus assumes a position such that the space charge of 
the ions in the accelerating region can reduce the applied electric field 
to zero. This space charge field is determined by the rate at which ions 
are transported across the plasma sheath. Since no field can be set up 
to reflect them, the ion current density is always equal to the maximum 

i 

available from the plasma. Since the plasma is (ideally) uniform, the 
ion current density is always uniform across the sheath and equal to the 
maximum available (assuming that sufficient potential is applied to the 
accelerator). Were it not for this property, the ion flow from the plas- 
ma sheath could be treated in identical fashion to the electron flow from 
a thermionic electron emitter. The principal property characterizing 
the plasma type ion emitter is thus its ability to adjust the boundary 
shape and position in order to satisfy the field and emission constraints 
discussed above. 

The geometry and potential profile on the axis of the unperturbed 
aperture studied are shown in Fig. 2. The emitting plasma density was 
assumed to be such as to give a nominal perveance of 3.4 x 10 ^ electron 
pervs. It can be seen that the ions, after extraction from the emitting 
plasma at a potential of approximately + 2.0 kV, are accelerated toward 
the accel electrode where the potential is - 2.0 kV . After passing through 
the aperture in the accel electrode, they are projected into an essentially 
field-free plasma at ground potential where neutralization takes place. 

The ion paths inside this latter plasma region are not deflected and the 
thrust from the aperture can be calculated from the ion trajectores at 
its (the downstream boundary) surface. 

3 . Downstream Plasma Boundary 

The boundary conditions applied for the solution of this 
downstream boundary are that the potential and the normal electric field 
be zero at the surface. The method used consisted of specifying the cor- 
rect potential of zero volts on a trial boundary and examining the normal 



0.090 in. 


E992-I0R2 



SCREEN ACCEL 

ELECTRODE ELECTRODE 



Fig. 2. Unperturbed electrode system and 
voltage profile on axis. 


component of electric field. The boundary was then deformed, until 
both these conditions were satisfied. This iterative method had to be 
employed since the potential and its derivative cannot be used as 
boundary conditions simultaneously as the Laplace and Poisson equa- 
tions to be solved would become "over specified . H 


The ion flow in a thruster consisting of a large number of closely 
packed co-planar accelerators can be considered to first order to be 
one dimensional. Hence, a first estimate of the location of the down- 
stream boundary was obtained assuming that the flow was between 
infinite plane grids as illustrated in Fig. 3. The electric field at the 
injection and exit plane of the ions was assumed to be zero since this 

is a boundary condition of plasma surface . In this case the relationship 

2 

between the spacings and voltages is given by 


_s 

d 


(1 - R -1 / 2 ) 


1/2 


(1 


+ 2R- 1 / 2 ) 


where 



and 

V = voltage of emitting plasma 

= voltage of neutralization plasma 

S = distance between emitting plasma and potential 
minimum 

d = distance between neutralization plasma and 
potential minimum. 


Hence a boundary at an axial distance of (d + S), as calculated by the 
above formulas, from the screen was used as the first estimate. This 
was successively refined until the normal electric field became negli- 
gible. The final solution of the downstream boundary was found to be 
located at an axial distance of 0.33 in. in front of the screen. 
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Fig. 3. Idealized system assuming planar symmetry. 
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It can be seen that the ion trajectories in the vicinity of the 
boundary are straight and hence the calculations of thrust are not 
sensitive to small movements of this boundary* Since changes in the 
upstream boundary have been found to produce negligible changes in 
the downstream boundary, the same downstream boundary was used 
for all solutions. 

4. Upstream Boundary 

The shape and position of the upstream plasma boundary 
exert a very strong effect on the ion flow and trajectories. In the 
calculations for the ion trajectories, the screen electrode was assumed 
to be at -40 V with respect to the emitting plasma and the accel elec- 
trode at - 4.0 kV with respect to the screen. In this case, the surface 
of the emitting plasma does not meet the screen electrode but will lie 
approximately parallel to it at a distance D (see Fig. 4). As a first 
estimate D was calculated by Child’s law for parallel diodes with a 
voltage difference of 40 V. 

The boundary conditions used for the solution of the emitting sur- 
face were the following: 

1 . The emitting surface must be an equipotential 

2. The normal derivative of the voltage must be 
zero across the surface 

3. The emission density at each point across the 
surface must equal the ion arrival rate at the 
plasma surface. 

These three conditions cannot be applied simultaneously in the 
solution of the Poisson equation and the trajectory equation without 
first determining the shape and location of the boundary. Hence an 
iterative method was used. 
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Fig. 4. Geometry at the emitter, showing the shield- 
ing effect of the screen. 
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The first of these conditions was met in the computer calcula- 
tions by assuming that the emitting surface is a Dirichlet boundary, 
and the second by assuming that the emission is space-charge limited. 
Since the ion current is determined solely by the ion arrival rate, the 
location and shape of the surface must be deformed until the space- 
charge limited emission is exactly maintained by the arriving ions. 
Hence the ion arrival rate across the surface must be calculated first 
from the desired perveance of the unperturbed beam. In this case the 
assumed perveance was 3.4 pperv.‘ As a result of the finite thickness 
of the screen electrode the arriving ion current density is not uniform 
across the screen aperture. For a plasma density of 10 ions/cm and 
an ion temperature of 300°C, the mean free path is long compared with 
the screen thickness. In this case the ion paths can be assumed to be 
collision free in the vicinity of the screen electrode and the ion arrival 
rate at a point on the emitter surface can be calculated from the angle 

0 shown in Fig. 4, which shows a diametral plane through the aperture, 
th. 

Hence for the i segment in this plane the current will be proportional 
to the angle 0^. The comparable angle from which the ions can arrive 
at the plasma boundary in the plane normal to that shown in Fig. 4 is 
4» . For the unperturbed case the current density in an element 
represented by one trajectory was calculated from the required perve- 
ance P and the product of the angles 0^ and 4^. using oc iji.fh and 



P N 

^ A s 


where N is the number of trajectories. 

Because of the inherent instability in the immediate vicinity of 
the emitter, a fundamental limitation arises in the accuracy to which 
the boundary can be located self- consistently {see Appendix I). As a 


The computer program is set up to calculate the perveance in elec- 
tron pervs* Hence this figure is given in these units. 
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result, the following method was used. A number of cases for vari- 
ous boundaries were first run as described above until the position of 
the boundary was known within approximately 0.0005 in. The best 
location for the boundary was then estimated and the correct current 
distribution is imposed upon it. This step is equivalent to assuming a 
small initial velocity of ions at the emitter, which is true in practice. 

The self-consistent plasma boundaries and beam profile for the unper- 
turbed geometry with a total beam current of 0.132 |jlA are shown in 
Fig. 5. 

5 . Calculation of Thrust for Unperturbed System 

The computer program, when used to solve problems 
with axial symmetry, divides the flow into a number of annular seg- 
ments centered on the axis of symmetry and calculates a self-consistent 
ion trajectory for each segment. Each trajectory so calculated illustrates 

the flow of all the ions emitted from its annular segment. If I. is the 

th 1 

current associated with the i annular segment, the total axial thrust 

T^ from the aperture is 

T . = kZI. n/V cos 6 . (1) 

A i i 

where the summation is for all the trajectories, k is a constant, and 

V = voltage difference between upstream and down- 
stream plasma 

4- = angle between trajectory and axis at downstream 

boundary (exit angle). 

The total axial thrust for the unperturbed geometry with a 
perveance of 3 . 4 x 10"^ was calculated to be 2.503 x 10 ^ lb. 

The transverse thrust T x is given by an expression similar 

to (1) 


T, = kZI.N/V sin i. ; 
1 i i 


( 2 ) 


because of axial symmetry, however, the summation leads to zero. 
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Fig. 5. Trajectories for unperturbed accelerator (case 1). 
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C . Thrust Calculations with Displaced Electrodes 

1 . Axial Displacements 

It has been found that the major contribution to the 
change in thrust for axial displacements results from the movement 
of the upstream plasma boundary rather than from the change of 
direction of the trajectories downstream. A change in screen-accel 
electrode spacing causes the form and position of the boundary to 
change, resulting in a change in emitting area and hence in the total 
emitted current. The total current I fc passing through a single 
aperture can be expressed as 

A 


i = y j (aa.) f. 

t l— J o 11 


(3) 


where J q is the density of current flowing through an element of area 
AA. on the emitting surface, F i is the "view factor" which accounts 
for the shielding effect of the screen electrode, A is the total curved 
surface area of the emitting plasma, and the summation is for all ele- 
ments on the emitting surface . It can be seen that an increase in area 

A will cause an increase in total current and hence an increase in the 
s 

thrust. The axial thrust from an aperture calculated from (1) can be 
expressed by means of (3) as 


T = k’s/v’ 2 j (AA.) F. cos <{>. . 

j O 11 1 


(4) 


From (4) we can find the effect of small changes in the various parameters 



3T. 

+ 2 7- 6(cos ^.) 

3(cos 40 1 


6(AA.) 6F. 6(cos </>■ ) 

= ST. A . - + ST. -=r— ■ + 2T. -r- 

i AA. i F. i cos p. 

li l 


(5) 
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The first two components result solely from changes in the emitting 
plasma boundary, while the third is an ion optical effect. The relative 
importance of the three components can be estimated from the magni- 
tudes of 


6(AA.) 6F. 6(cos 4 .) 

1 1 , 1 

Aa. 9 F. an cos 

ii l 


When the solution obtained for the unperturbed system is compared with 
that for a 22% decrease in accel- screen spacing, the change in thrust 
resulting from ion optics can be shown to be small compared with that 
resulting from the boundary movement. This can readily be seen if 
6T^ is expressed as 


6T a 



6(cos </> 


b 


= ZT. 


l 



+ ZT. 

l 


6(cos ^.) 
cos d. 

l 


( 6 ) 


and if the maximum values 6I./L and 6(cos >(.)/cos are compared. 
Since 

(cos /.) 

6 t — = tan 6«(. , 

cos p. 1 1 

i 

is small (less than 20° for all trajectories and tan 20° = 0.36), and 
the maximum value of 6^ < 0.1, this term will be less than 0.036. 
However, the maximum value of 6I./L = 0.13. Hence the change in 
thrust resulting from the change in current is the dominant factor in 
this case; this result is somewhat unexpected at first glance. 
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When the actual values rather than the maximum values are used 
for the 22% decrease in axial spacing, summing for all the trajectories 
we obtain 


hence, for 6T^, 


6L 


= 0.117; 


6L , 

ST. -s-— = 2. 9 x 10“ lb 
1 I. 

1 


6(cos <£.) _ 

ST. r— = 3.6x10' lb 


1 cos <f>. 

1 


However, the increase in emitted current results from two factors: the 
increase in emission area and the increase in view factor. We can 
therefore write 


I, = SI. = J 0 2(A Ai ) F 
6I t = 2 WXJ 6(AA i» + S 7F 6F 

= ZJ 6(AA.) F + SJ 6F (AA.) 

O 1 O 1 

6(AA.) 

= SJ o (AA.) F + SJ o 6F (AA.) . (7) 


Since the trajectories calculated represent annuli of equal thick- 
ness, the fractional change in emitting area will be the same for all ele- 
ments and we can write 


6(AA.) 

aa. 


(is constant) 


6A 

s 

A 

s 


( 8 ) 
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Therefore, 


6A 

61. = I. — + 2J 6F.(AA.) 
t t A o i* r 

s 


or 



6A SJ 6F.(AA.) 
s ^ o i x V 

A L 

s t 


( 9 ) 


In this case 6A /A was calculated to be approximately 0.09; hence 
s s 

the change in current resulting from the change in view factor is 
approximately 2.7%. 

We can now write the change in thrust 6T. as 

A 






( 10 ) 


where 

6Tj A - change in axial thrust due to change in ion 
paths alone 

6T^ A = change in axial thrust due to change in emitting 
area 

6T.j A = change in axial thrust due to change in screen 
electrode shielding. 

A change in the ion arrival rate will also introduce a perturbation which 
preserves the axial symmetry of the system. In this case the change in 
thrust is given by 
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( 11 ) 


6T a 


ST. 


61. 

1 


1 I. 
1 


+ ST. 
i 


6(cos ^.) 
cos p. 


where 

6A 6j 

61. = I. -r- 5 - + I. — + 2 j 6F (Aa. ) (12) 

i t A t — J o r 

s J 

J o 

because the ion arrival rate j varies. 

J o 

Using the above notation we can write 







(13) 


where 6T . A = the change in thrust because of the change in ion 
4 , A 

arrival rate . 


Five cases with axial symmetry have been solved for three 

"ion arrival rates, M and the total thrust is given in Tables I, II, and III. 

The three values of ion arrival rate assumed in Table I are such as to 

maintain an ion current density at the screen aperture of the unperturbed 

Z Z 

system of 7.87 mA/cm (cases 1 and 2), 4.5 mA/cm (case 3) and 
2.28 mA/cm (cases 4 and 5). The electrode perturbation considered 
was an axial movement of the accel electrode of 0.020 in toward the 
screen. The trajectories calculated for cases 1 through 5 are shown 
in Fig . 5 to 9 . 

For small perturbations the change in thrust can be divided into 
components parts according to eq . (13). This has been done for the 
cases 1, 2, 4, and 5 shown in Table I, and the results are shown in 
Table II. 
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Fig. 6. Ion trajectories for case 2. 
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Fig. 7. Ion trajectories for case 3. 
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ACCEL ELECTRODE 


Fig. 8. Ion trajectories for case 4. 
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Fig. 9. Ion trajectories for case 5. 
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TABLE I 


Thrust Calculated for Electrode Configurations with 
Axial Symmetry in Unperturbed System 


Case 

Number 

Screen/ 

Accel 

Spacing, in. 

Average 
Current 
Density^ 
mA/ cm 

Electron 
Perv 
X 10" 6 

Axial Thrust, 
lb x 10“ 5 

Constant Ion Arrival Rate = 3.4 

x 10’ 6 electron perv 

i 

0.090 

7.87 

3.4 

2.503 

2 

0.070 

8.87 

3.83 

2.829 

Constant Ion Arrival Rate = 1 . 942 x 10 ^ electron perv 

3 

0.090 

4.5 

1 . 942 

1 . 502 

Constant Ion Arrival Rate = 0. 985 x 10 ^ electron perv 

4 

0.090 

2.28 

0.985 

0.821 

5 

0.070 

2.56 

1.105 

0. 887 

2^ 

Evaluated at screen aperture. 


TABLE II 


Change in Thrust due to 0.020 in. Movement of Accel 
Electrode in the Unperturbed System 


Case 
Numbe r 

Axial 
Thrust,, 
lb x 10- 5 

6T/T, % 

6^/T, % 

6T 2 /T, % 

6T 3 /T, % 

Ion Arrival Rate = 3 . 4 x 1 0 electron perv 

l a 

2.503 

— 

— 

— 

— 

2 

2.829 

13 

1.43 

8.9 

2.67 

Ion Arrival Rate = 0. 985 x 10"^ electron perv 

4 a 

0.821 

— 

— 

— 

— 

5 

0.887 

8.0 

3.4 

4.5 

0.1 

Considered as the unperturbed system. 


The above analysis applies strictly to small perturbations only, 
since it assumes linearity neglecting products of first order terms. An 
estimate of the changes in thrust due to large reductions in ion arrival 
rate in the unperturbed electrode system can be obtained if we write 

6T a = K'Jv |si. 2 cos (d i2 ) - SI U cos (s^) 

= K*Jv js(I n + 61.) cos (J. + 6^.) - 21. j cos (j^) 

and 


6T A1 = K'Vv 2 I u cos fczf.j + Sj^) - cos (^.j) 
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This assumes 261 cos + 6^) = 61 cos ^ 


6T 


9T 


A2 ~ & (A A ' )' 6{AA s } 


9t 

6T A3 = W7 5F i 


9T o 

-A4 ' ^ *, • 


The results broken down in this way for the three ion arrival rates are 
shown in Table III. (It can be seen that the error in this assumption is 
about 4%) . 


TABLE HI 

Thrust Calculated for Accelerator with Axial Symmetry 


Case 
Numbe r 

Axial 
Thrust, 
lb x 10- 5 

6T/T, 
’ % 

6T i/T, 
% 

STt/T, 

% 

6T 3 /T, 

% 

6T 4 /T, 

% 

1 

2. 503 a 

— 

— 

— 

— 

— 

2 

1.502 

-40 

2.9 

8.6 

2.7 

- 54.3 

3 

0.821 

! -67.2 

5.2 

14. 

9.5 

-71.3 

2L 

Considered as the unperturbed system. 


From Fig. 10 it can be seen that the variation in thrust with 
emitted current is linear and hence 
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NOMINAL PERVEANCE IN UNPERTURBED SYSTEM, xIO"® ELECTRON PERVS 


Fig. 10 


Change in thrust due to ion arrival rate with constant screen 
to accel spacing of 0. 090 in. 




9T 

a T 


change in thrust 
change in ion arrival 


can be evaluated as 


= 0.01725 mlb/mA 


similarly, the derivative 


dT (%) 
dd 


percentage change in axial thrust 
percentage change in electrode spacing 


at constant 
ion arrival 
rate 


3T. 

= 2 -ffizskj) 


9 T. 

+ s 


+ 2 


9T i 

Tp: 

1 


can be evaluated for two ion arrival rates. 

From Fig. 11 and 12 we see that 

0.103, %/mil 

atP = 3.4x10 
electron pervs 


dT<%) 

dd 


dT (%) 
dd 


= 0.419 %/mil . 


at P s 0.985 x 10 
electron pervs 


From the above results it is possible to estimate the thrust 
for small perturbations which preserve the axial symmetry of the system. 
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THRUST, lb* 10 



PERCENTAGE OF UNPERTURBED THRUST 











2 , 


Transverse Displacements 


When the accel-electrode is displaced transversely, 
the axis of symmetry of the aperture is destroyed; hence the ion 
motion cannot be solved assuming circular symmetry. 

A transverse misalignment is illustrated in Fig. 13, where the 
axes of the screen electrode and the accel electrode are shown dis- 
olaced by a distance ea (where a is the radius of the accel aperture). 

The R - Q equation of the displaced surface can be found by 
expressing 


R = x + iy 


(14) 


where 


x = acos 0 
y = ea + a sin 0 ; 


hence 


R = (x 


, 2 . 1/2 
+ y ) 


(15) 


or to first order 

R = (a 2 + 2ea 2 sin O) 1 ^ 2 = a {1 + e sin 0) . (16) 

It can thus be seen that this type of displacement introduces a per- 
turbation in the radial coordinate R which varies as sin 0. 

If two sections are drawn through the system normal to each 
other, representing the planes of least and greatest perturbations,, the 
axis lines of the two electrodes will coincide in the first plane and be 
displaced by 6 a in the second plane. The trajectory equations can 
be solved as planar problems for these two cases and the results will 
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be as illustrated in Fig. 14. It can be seen that the calculated trajec- 
tories will differ. The effect of this transverse displacement can be 
thought of as displacing each trajectory through an angle 6. . Hence, 
for the trajectories nearest the screen the change in exit angle is 6 ^ 
at the top and 6^ at the bottom. 

If it is assumed that the change in exit angle varies as sin 0 in 
the azimuthal direction (as shown in Fig. 15), the angular displacement 
of the outermost trajectory in a plane at an angle (tt/ 2 - 0) from the 
plane of the displacement can be written as 


" 7 (6 1 + 6 27 ) sin 0 + I ( 6 i - 6 27 ) * (17) 

Hence in a plane normal to the displacement where 0 = 0, the change 
in exit angle 6?T is l/2 (6^ - ^27^* P^ ane displacement 

at 0 = tt/2, b</> ^ at 0 = - tt/2, 6^ = 6 ^ . Hence from a small 

element (shown in Fig. 20) where the current density is j ^ the axial 
thrust is given by 


AT. = j. 6S. cos t. , s/V (18) 

i,A J i l i i v ' 

where 1 is the exit angle 6S. = R. 60 

l l i 

R. = radial coordinate of the element 

l 

60 = angular width of the element. 


The refore , 


AT. a = R.t. j . cos 60 . 

l, A ii J i l 

Hence the total thrust from the annulus shown is 


(19) 


37 


E992 *12 



Fig. 14. Effect of translational misalignment of the ion paths calcu 
lated assuming planar symmetry (dotted lines represent 
perturbed solutions). 
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Fig. 15. Three-dimensional representation of the perturbed and un- 
perturbed trajectories. 
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T. A = 2(R.t.j. ) n/v cos 60 . 

l, A 1 i J i i 


( 20 ) 


In this equation we can write J^O = Rjt.jj 60. J\ is the current per 
radian in the i^ trajectory; therefore, J\6Q is the current contained 
in an element of the i th trajectory subtended by an angle 60 at the 
axis. Hence I. = 2J i 60. Since boundary conditions for the upstream 
plasma are that the emission density be equivalent to the ion arrival 
rate at the surface and that the ion density be assumed constant across 
a single aperture, it can be shown that variations in J. with 0 will 
produce at most second order perturbations in the thrust. (Variations 
in J. will be produced by variations in the shielding effect of the 
screen electrode only because of the asymmetry of the upstream 
boundary.) Hence J = (L/2 tt). In (20), 

</> ' = -4 (6. + 1 >. ) sin 0 +4 (6. - 6.) (21) 

’l *i 2 l l 2 l i 

whe re 


^ . = exit angles for the unperturbed axisymmetric run 

5. = maximum positive angular displacement 

"6. = maximum negative angular displacement. 


Hence 


T. 

l 


I.'v/V 

l 

2tt 



( 6 . + 6 .) 

l l n 

+ ^ sin 9 


d0 


I.'s/V 






( 6 . + 6 X ) 

2 


sin 6 


dO 





( 22 ) 
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However, 



sin 0 


dO = 2ir J 

o 



+ 6 . 

1 

2 


■) 


where J is the Bessel function of order zero and 
o 


2 J / 6 . + 6 . \ 

/ sin l 2 sinOy 


d0 = 0 . 


Hence 


/6. + 6.\ f 6. - 6. \ 


The total perturbed axial thrust is therefore given by 


/ 6 . + 6 . \ / 6 .- 6 .\ 

T o = 4 ^ J o (rirv c ° s ('i + -4-0 


6 . - 6 . 


(23) 


(24) 


The self-consistent ion trajectories solved assuming a 0.0142 in. 
transverse displacement and planar symmetry and a normal ion arrival 
rate equivalent to a perveance of 3.4 jipervs are shown in Fig. 16. A 
self-consistent upstream boundary was found for each case. It can be 
seen that the beam is deflected in an opposite direction to the perturba- 
tion, and that the two beam edges are moved in this direction by 
approximately equal amounts. This result is in agreement with pre- 
vious work. The angular deflection at the edges is about 4°; hence 
for the beam edge the angular deflection can be calculated to be 
0.028°/mil. However, the mean angular deflection *6^ is calculated 
as 



( 6 ^) 

“ N~ 
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Fig. 16. Ion trajectories calculated assuming planar symmetry for 
case 6. 
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whe re 


N 


number of trajectories 

8 . 2 ° . 


Hence if it is assumed that the mean angular deflection is a linear 
function of transverse displacement of the accel electrode, the pro- 
portionality constant can be calculated as 0.578°/mil. In metric units 

this is equivalent to 23. 2°/mm, which compares with a value of 29°/mm 

4 

reported by other workers. This difference can be explained by a num- 
ber of factors. The geometry and voltages assumed in this case are 
somewhat different. Thus the lens produced in our geometry is 
slightly stronger (which tends to reduce the divergence) and in previous 
work no account was taken of the relocation of the emitting plasma 
caused by the electrode perturbation; in addition, the correct position 
of the downstream boundary was not determined. 

If the total thrust is calculated using for this case 


T = SI. 

l 


\Fv j 

o 



6 . 




a value of 2.47 2 x 10 


lb is obtained. 


The transverse thrust T' is a vector which can be represented 
by its components and T!> where T'^ is the component of thrust 

parallel to the displacement €a and is the component normal to it. 

Then 

I. 

*v/v sin ^ sin OdO (25) 



and 



I. 

\Tv sin ^!! cos OdO 


(26) 
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where fT = ^. + (l/2) (6. + 6.) sin 0 + (l/2) (6. - 6.) and the summa- 
tion is for i = 1 to 27 . 

We can write 


T'j = 2 


2 ? I. _ ( / 6.-6! / 6. +6. \ 

J 2^ ^ j sin [fi + ~~~2 cos ( 1 2 1 sin0 J 


/ 6 .- 6 . / 6 . + 6 . 

+ cos ( + — x - 2 — — sin( 1 2 — — sin 0 ) } sin OdO (27) 


and 


TL 


y 

} I. ( / 6. - "Fa /6. + 6. \ 

= Z J i + 1 2 J cos ( - 1 2 sinQ ) 


/ 6 . - 6 . \ / 6 . + 6 . ^ 

+ cos 2 — — ) sinf — ^2 — “ sin ojl) cos OdO 


(28) 


Noting that 
2 tt 


J sin ("Z ^ 6 i + V sin y 


cos 0=0 


2 tt 



i _ 

\ f b 


/ sin (i 

> (6 i + 6^) sin Q 

J sin 6 = 2n Jj ( — 

4-0 


2tt 

j cos^| (6^ + 6.) sin 0^ cos 0 = 0 


o 

2tt 


J cos (6^ + 6.) sin 0^ si 


sin 0=0 , 


(2 9) 
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we can write 


21. *s/V 

L 





( 30 ) 


Ti, = 0. (31) 

The values of thrust have been calculated for both a 0.0142 in. 
and a 0.0071 in. transverse displacement of the accel electrode for 
two values of ion arrival rate. The ion trajectories calculated are 
showninFigs. 16 to 19. The results are given in Table IV and are 
illustrated in Figs. 20, 21, and 22. 

TABLE IV 

Thrust Calculated for Accelerator with Transverse Perturbations 


Case 
Numbe r 

Nominal 

Perveance, 

Electron 

Pervs 

T ransverse 
Perturbation 
of Accel 
x 0.001 in. 

T?2 

Average 

Deflection, 

deg 

Axial 

Thrust, lb 

Transverse 
Thrust, lb 

Rotation 
of Thrust 
Vector, deg 

6 

3.4 x 10* 6 

0.0071 

4. 13 

2.487 x 10* 5 

1.278 x 10' 6 

2.94 

7 

3.4 x 10" 6 

0.0142 

8.2 

2.472 x 10' 5 

2.529 x 10' 6 

5.84 

8 

0.985 x 10' 6 

0.0071 

4.49 

8. 12 x 10' 6 

3.573 x 10' 7 

2.52 

9 

0.985 x 10' 6 

0.0142 

8.8 

8.032 x 10' 6 

7.085 x 10* 7 

5.04 


It can be seen that for the ranges of displacement considered the change 
in thrust is a linear function of the electrode spacing. 

3. Tilt 

For an engine with a diameter of 30 cm and a spacing 
between the screen and accel electrode of 0.090 in., the maximum 
angle for gross tilt is of the order of tan ^ 0.007. This is equivalent 
to an angle of about 0.4°. (However, localized warping could increase 
the angle to some extent.) A case has been run with a symmetric tilt 
of 1.5° about the unperturbed center line. This produces a maximum 
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Fig. 17. 


Ion trajectories calculated assuming planar symmetry for 
case 7. 
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Fig. 18. Ion trajectories calculated assuming planar symmetry for 
case 8. 


47 





Fig. 19. Ion trajectories calculated assuming planar symmetry for 
case 9. 
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TRANSVERSE THRU SI 



DISPLACEMENT, in. 

Fig. 22. Transverse thrust due to displacement of accel electrode 




displacement of the a c cel- electrode of 0.002 in. The results from this 

run show no significant change in ion optics or emission density. This 

4 

is in general agreement with the other reported work, where it was 
found that for tilt angles of less than 5° no meaningful results could be 
obtained since the small calculated beam deflection approaches the 
accuracy of the computer program. 

However, tilt may be a significant factor at the edge of the 
engine where the displacement of the accel-electrode can be regarded 
as a combination of an axial movement and a tilt. In this case a satis- 
factory solution can be found by considering the appropriate axial dis- 
placement alone. 

D. Conclusions 

The iqn optics have been solved self- consistently for ten differ- 
ent configurations and' the total thrust from a single aperture calculated 
for each. The thrust from the unperturbed system with a nominal per- 
veance of 3.4 x 10~ electron pervs was calculated to be 2.503 x 10 lb. 
Three types of perturbation have been considered: (1) changes in ion 
emission rate, (2) changes in the axial spacing of the electrodes, and 
(3) transverse misalignments. The first two types can be solved directly 
since they both possess axial symmetry. The third type produces a dis- 
turbed flow about a curvilinear axis and cannot be solved precisely by 
presently known techniques. A first order solution has been formed for 
this . 


It has been found that changes in thrust due to perturbations 
which preserve the axial symmetry result primarily from the reloca- 
tion of the emitting plasma, and the influence of ion optics has been 
found to be small. However, for the transverse misalignments con- 
sidered (0.0071 in. and 0.0142 in.), the effect of the upstream boundary 
movement is small and the change in ion optics is more significant. 

The results can be expressed in terms- of the parameters defined 

below. 
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a 


P 


change m axial thrust 

hange in axial spacing of electrodes 


at constant ion emis- 
sion rate equivalent to 
an electron perve- 
ance of P 


change in axial thrust 
change in emitted current 


v 



change in axial thrust 

transverse displacement of accel 


at constant ion emission 
rate equivalent to an elec- 
tron perveance of P 


*1 



change in transverse thrust 
transverse displacement of accel 


at constant ion emission 
rate equivalent to an elec- 
tron perveance of P 


at constant ion emission 
rate equivalent to an elec- 
tron perveance of P 

The values derived for the above constants are given in Table V. These 
values will be used in Section IV, together with the calculated magnitudes 
of electrode displacement, to provide the estimated total thrust vector 
variation of a single thruster. 


co 


rotation oi tnrust vector 

transverse displacement of the accel 


TABLE V 

Parameters Defining Change in Thrust Vector for Various Perturbations 


P (Nominal 
Perveance), 
Electron 
Pervs 

lb/0. 001 m. 

lb/mA 

Y | p 

lb/0.001 in. 

^Ip 

lb/0.001 in. 

Hp 

deg/0.001 m. 

V 

3. 4 x 10" b 

- 1 63 x 1 0 7 

1.765 x 10" 5 

Z. 18 x lo" 8 

-1 78 x 10" 7 

0.411 

0.985 x 10" 6 

-0.33 x 10" 7 

0.684 x 10" 5 

1.26 x 10" 8 

-0.498 x 10" 7 

0.354 
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in. CAUSES AND MAGNITUDES OF ELECTRODE MISALIGNMENT 

A. Introduction 

Under ideal conditions the geometric composition of a Kaufman 
thruster may be described as an arrangement of right circular cylinders. 
The main thrust axis is defined by the large cylindrical discharge chamber 
with the circular electrodes lying in parallel planes at one end. Thousands 
of hole pairs in the electrodes, located symmetrically about the thrust 
axis, form an array of small cylindrical accelerating structures whose 
axes are. all aligned parallel to the thrust axis. The over-all thrust 
vector produced passes precisely through the geometric center of the 
electrodes and is parallel to the common sixes of all the cylinders; no 
torques are produced about any axis. An additional assumption required 
for this ideal thruster operation is that the ion density in the discharge 
chamber be axisymmetric. If this assumption is violated or if any of 
the aligned cylindrical geometries become disturbed, perturbations in 
the thrust vector will be incurred. 

Of the many possible causes of thrust vector perturbation it is 
reasonable to study those which are most likely to occur and which pro- 
duce the largest thrust disturbances. The thrust vector disturbances 
produced are analyzed so that tolerance can be specified for the current 
design efforts. 

The cause of thrust vector perturbation which was analyzed in 
sufficient detail for quantitative prediction is the thermal buckling of 
the electrodes. In contrast to many other causes of thrust vector dis- 
turbance, the misalignments in the grid system which result from heat- 
ing (an assembly initially aligned at room temperature) to operating 
temperatures are inherent in the basic thruster design. The largest 
and most often observed thermal electrode motion is axial displacement 
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(bowing). Because of the frequent occurrence and seriousness of the 
buckling problem in operating thrusters, electrodes are usually given 
an initial bow so that the direction of subsequent motion will be deter- 
mined and preset ,to coincide. One may suspect that axial displace- 
ments equal to half the initial grid spacing are common. It will be 
shown in the subsequent sections that a spacing change from a uniform 
0. 090 in. to a parabolic distribution, which varies from 0. 090 in. at 
the edge to 0. 045 in. at the center, will cause an 18% change in thrust 
magnitude. 

Although electrode spacing changes will be the only cause of 
misalignment that will be modeled analytically, the thrust vector 
perturbations caused by a number of modes of electrode misalignment 
are evaluated for thrust vector disturbance. 

The thermal buckling analysis is presented in its two main parts: 
(1) the thermal analysis leading to the grid system temperature distri- 
butions, and (2) the stress analysis of the mechanically loaded electrodes 
with given temperature distribution. 

These analyses will yield the expected electrode perturbations, 
which will be combined (in Section IV) with the coefficients, relating 
the thrust vector to electrode displacement, to give the expected thrust 
vector variation. 

B. Thermal Simulation of the Grid System 

In order to provide a guide to initial thrust analysis, before a 
thermal simulation of the grid system was attempted, temperature 
distributions were measured on a thermal mockup of a bombardment 
thruster. A schematic drawing of the mockup is shown in Fig. 23. 

Each spot represents the location of one of the 24 thermocouples. The 
heat sources were simulated by heating coils located at the front of 
anode A, over the remaining anode surface B, and at the cathode 
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position. Three temperature distributions were measured for various 
combinations of input heater powers, with the mockup in a vacuum tank. 
TypicaL data are also shown in Fig. 23. 

Efforts to predict these data by using analytical models based on 
nonreflective radiative transfer were unsuccessful. Upon further con- 
sideration of the data (particularly the rather uniform equilibrium 
temperature distributions across the electrodes), it was suspected that 
a large number of multiple reflections occur inside the discharge 
chamber; accordingly, the heat flux distribution across the screen 
electrode is fairly uniform. Because it was recognized that a difficult 
multiple reflection radiation problem was at hand, the simulation was 
directed toward a detailed analysis of the grid system only, under the 
assumption of various incident heat flux distributions from the discharge 
chamber. 

The multiple reflection heat transfer in the grid system was 
simulated by the existing thermal analysis computer program TAS-1B. 
This is a digital computer program that solves for the steady- state 
temperature distribution in a lumped parameter network of temperature 
points (nodes) and heat flow paths (resistors). A node is provided for 
each constant temperature portion of the device being modeled. The 
network may contain up to 80 nodes and each node may be connected 
to any or all of the others. Any two nodes may be connected by two 
resistors. One represents conduction and/or convection; the other 
represents infrared radiation. In addition to interchanging heat with 
other nodes, each node may receive radiation and/or a direct heat input 
(such as from internal dissipation of electrical energy). 

Considerable effort was devoted to formulating a thermal network 
for the two perforated electrodes. The definition of nodes and the sub- 
sequent calculation of their view factors and thermal conductivities is 
discussed below in Section III-B-1. In Section III-B-2 a parametric 
study is given of electrode temperature distributions which are calcu- 
lated from the thermal simulation. 
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1 . 


Network Model 


A thermal model, of the grid system was developed in the 
form of a network of nodes, each of which corresponds to a constant 
temperature region on an electrode. For axisymmetric heating the 

t 

nodes N^ and are washer shaped (annuli) regions on the electrodes 
(see Fig. 24). Heat is transferred between adjacent washers of each 
electrode by conduction and between opposite washers of the two 
electrodes by radiation. An additional constant temperature node Ng, 
which represents the background, absorbs heat radiated directly from 
the accelerator electrode and from the screen electrode by way of the 
holes in the accelerator electrode. The heat input to the grid system 
from the discharge chamber is represented by a distribution of black 
body radiation which is incident on screen electrode and varies radially 
in intensity. The thermal network corresponding to this model is 
shown in Fig. 25. The node Nb is fixed at a constant temperature 
(0°F in our examples) and is assumed to be a black body (e = 1) whose 
area is much larger than any other node. Under these assumptions 
Nt, may be regarded as a heat sink. For given values of incident heat, 

j 

flux Q^, the temperatures of the other nodes N^, can be calculated 
once the node coupling by the conductance paths C(i, j) and radiation 
paths R(i, j) is established. The evaluation of the thermal conductivities 
between the nodes of the model is straightforward and will be considered 
first. The evaluation of the geometric (view) factors which govern the 
radiative heat transfer in the model is more complicated; only a brief 
discussion will be given which is based on the detailed view factor 
calculations in Appendix I. 

a. Heat Conductance in the Electrodes — The tempera- 

ture of a (washer shaped) node on an electrode shown in Fig. 26(a) will be 
taken as the value of the temperature at the circle halfway between the 
edges of the washer. Thus, the temperatures of nodes N^, N^, and 
are taken at the locations R^, R^, and R^. The thermal conductance 
between nodes and + 1 is then equal to the conductance between the 
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inside and outside edges of the washer shown in Fig. 26(b). The formula 
for this conductance is 


C. 


i+1 


ZirtK 1 


log 


R i+1 

R. 

1 


where 

t = thickness of the washer 

K* s thermal conductivity of the washer. 

For a perforated plate of open area fraction f and material thermal con- 
ductivity K, we ’take K 1 = (1 - f) K. 


b. View Factors for Radiative Heat Transfer in the 

\ — - • — — — 

Grid System — Because of the small spacing 
between the electrodes relative to the radii of the annular nodes (see 
Fig. 24), it is reasonable to assume that radiative heat transfer occurs 
only between the directly opposing nodes (i. e. , between N. and 1NL). 
Under this assumption the desired view factors for the grid system can 
be calculated on the basis of two infinite parallel perforated plates. 

. i 

The geometry of a cross section of two aligned perforated plates S and 
A is shown in Fig. 27. The nodes for which view factors are to be 
calculated are indicated by the numbers 1 through 7. The actual cal- 
culations are carried out in Appendix II; because of their complexity, 
approximations are made. The nodes 1 through 6 describe in detail 

T 

the coupling between each set of nodes (N. , N.). The over-all detailed 
network for the grid system has three sets of six nodes each, as shown 
in Fig. 28. 

The heat incident on the screen electrode is first divided (into 

Q^, Q 2> and Q^) between the three nodes N^, N£, and A radial 

heat flux distribution is assumed which varies linearly between Q q at 

the center of the screen and Q„ at the edge. The heat components Q. 

-K l 

are then split into (1 - f g )Q.. U n P u t to the material surface of.N. (node 

3)) and f' Q. (input to the open surface area of N. (node 1)). Here f is 
si is 

the open area fraction of the screen. In the heat transfer calculations 
it is assumed that node 1 is a black body (e =1). 
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2 . 


Grid System Temperature Distributions 

A parametric study of the equilibrium grid system tempera- 
ture distributions was made for a variety of electrode spacings and input 
power profiles. The calculations were carried out by the computer 
program T.AS-1B and employ the network model shown in Fig. 28. The 
nodal areas and thermal conductivities, which correspond to the assumed 
electrode geometry given in Table VI, are given in Table VII. The view 
factors corresponding to three values of electrode spacing 0. 09 in., 0. 07 in., 
and 0. 05 in. are given in Table VIII. 


TABLE VI 

Assumed Electrode Geometry 



Screen 

Accelerator 

Diameter 

30 cm 

30 cm 

Thickness 

0. 030 in. 

0. 060 in. 

Open area 

70% 

40% 

Hole diameter 

0. 187 in. 

0. 142 in. 

Hole spacing 

0. 213 in. 

0. 213 in. 


The sensitivity of the electrode temperature distributions to 
variations in (a) electrode spacing, (b) discharge power distribution, 

(c) total discharge power, and {d) interception power was studied and 
will now be discussed. 

, a. Electrode Spacing — The sensitivity of the electrode 

temperatures to the spacing between the electrodes is due to the depen- 
dence of the view factors on spacing (see Table VIII). 

The six curves in Fig. 29 give the screen and accelerator 
electrode temperatures for three values of spacing. As shown, 
decreasing the interelectrode spacing lowers the over-all temperatures 
of both electrodes. Apparently, the grid system appears more trans- 
parent to the incident radiation as the electrode spacing decreases. 
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ACCELERATOR TEMPERATURE, °C SCREEN TEMPERATURE, °C 


E992 - 



DISTANCE FROM ELECTRODE CENTER, cm 


Fig* 29* Dependence of grid system temperature distributions on electrode spacing 


6 ~ 


Since the actual decrease in temperatures is small and the temperature 
gradients across the electrodes remain unchanged as the spacing varies 
the assumption that the grid system temperature distributions can be 
calculated independently from the thermal -warping is substantiated. 

TABLE VII 

Nodal Areas, Power Inputs, and Thermal Conductivities 


Areas 


= 8.55 
= 5 . 50 
= 3.55 
= 7 . 30 
= 8.25 
= 7 . 30 
= 10 7 
= 25 . 6 
= 16.4 


= 21 . 6 
= 24 . 6 
= 21 . 6 
= 43 . 0 
= 27 . 5 


= 36 . 2 
= 41 . 7 
= 36 . 2 


Thermal Conductivities, W/°C 


'3, 10 


' 10 , 16 


' 6 , 13 


' 13,19 


0 . 314 


0 . 674 


1 . 28 


2 . 74 


C 2 , 3 " C 9 , 10 " C 15 , 16 


C 4, 5 " C 17 , 18 " C 5, 6 


C 12 , 13 " C 18 , 19 

..6 W 


a Instead of calculating the very high conductances be- 
tween these nodes, we insured that they arrive at equal 
temperatures by coupling them with very high con- 
ductances. 





TABLE VIII 

Grid System View Factors 



Electrode Spacing 


0. 90 in. 

0. 070 in. 

0. 050 in. 

F l, 2 ~ F 8, 9 ~ F 14, 15 

0. 27 

0. 27 

0. 27 

F l, 4 = F 8, 11 = F 14, 17 

0. 38 

0. 35 

0. 33 

1,5 * 8, 12 r 14, 18 

0. 14 

0. 16 

0. 17 

F 1, 7 = F 8,7 = F 14, 7 

0. 21 

0. 22 

0.23 

F 2,2 ■ F 9,9 = F 15, 15 

0. 15 

0. 15 

0. 15 

F 2,4 = F 9, 11 = F 15, 17 

0. 31 

0.29 

0. 26 

F 2, 5 = F 9, 12 = F 1 5, 18 

0. 05 

0. 07 

0. 09 

F 2, 7 = F 9, 7 F 15, 7 

0. 05 

0. 06 

0. 08 

F 3, 4 = F 10, 11 = F 16, 17 

0. 60 

0. 68 

0. 75 

3,5 £ 10, 12 r 16, 18 

0. 3 

0. 26 

0. 21 

F = F - F 

J 3, 7 r 10, 7 £ 16, 7 

0. 1 

0. 06 

0. 04 

F 5, 5 = F 12, 12 = F 18, 18 

0. 34 

0. 34 

0. 34 

F 5, 7 = F 12, 7 = F 18, 17 

0. 33 

0. 33 

0. 33 

F 6, 7 = F 13, 7 = F 14, 7 

1. 0 

1.0 

1. 0 











b. Discharge Power Distribution — As shown in Fig. 24, 
the heat flux distribution from the discharge chamber is assumed to vary 
linearly between the values at the center and at the edge. Since heat 
conduction is small in the screen electrode (because it is thin and highly 
perforated) its temperature gradient, from center to edge, will depend 
strongly on the ratio Q^/Q.^. This dependence is illustrated by the three 
sets of temperature distributions in Fig. 30, which correspond to Q /Q t> = 1, 
1. 5, 2. For all three cases the total discharge power incident on the screen 
is 250 W and the interception power into the accelerator is 20 W (iO W into 
node N, and 10 W into N_). The heat intensity ratio Q /Q D = 1.5 gives rise 
to a 40°C temperature variation from the screen center to edge. This 
value is considered reasonable and is used in the subsequent sensitivity 
studies. 

c. Total Discharge Power — The total heat flux incident 
on the screen electrode from the discharge chamber is highly dependent 

on the specific discharge chamber performance characteristics (e. g. , the 
ionization efficiency) and the operating point of the thruster. The variations 
in the screen and accelerator center temperatures and center to edge 
temperature differences have therefore been calculated as functions of 
the total discharge power. The curves are shown in Fig. 31. For a 
total discharge chamber power (cathode power plus arc power) of 400 W 
it is reasonable to expect that half (200 W) will be incident on the screen. 

In this case the electrode temperatures are the following: 

Screen Accelerator 

Center Temperature 430°C 290°C 

Edge Temperature 407°C 269°C 

d. Interception Power — As illustrated in Figs. 30 
and 31, the effects of varying the discharge power intensity ratio Q /Q 

O XV 

and the total discharge power are most noticeable on the screen electrode. 
Another source of power into the grid system is the heating of the 
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accelerator by intercepted ions, which may be produced by charge ex- 
change collisions. Relatively small amounts (10 to 2.0 1AT) of interception 
power cause noticeable accelerator temperature changes, as shown in 
Fig. 32. The interception power is divided equally between nodes Nj 
and n! (see Fig. 32). Since the areas of nodes N., N.,, N 3 are in the 
ratios 1:4:9, the interception power distribution assumed is peaked in 
the center. The accelerator center to edge temperature difference is 
increased from 10°C to 30° C by 20 W of interception power. 

In Fig. 33, a breakdown of the heat flow paths in the grid system 
is shown for the case of 200 W incident on the screen in the distribution 
Q /Q r = 1.5 and 20 W of interception power. The quantities of heat 
indicated represent the net heat transferred after multiple reflections. 

Of the 29. 2 W incident on the node Nj, 20. 8 W goes to space (part by 
radiation and part by transfer straight through the holes), 6. 4 W is 
absorbed by the accelerator node N^and 2 W are conducted to the 
adjacent screen node N,. 

3. Correlation of the Thermal Model with Experiment 

An opportunity to correlate the thermal model which we 
have developed with some experimental data was provided by some 
recent JPL temperature measurements on an operating thruster. The 
electrode geometry of the. J PL thruster is the following: 



Screen 

Accelerator 

Diameter 

20 cm 

20 cm 

Thickness 

0. 030 in. 

0. 110 in. 

Open Area 

70% 

40% 

Hole Diameter 

0. 187 in. 

0. 142 in. 

Hole Spacing 

0. 213 in. 

0. 213 in. 


and the electrode spacing is 0. 070 in. The corresponding nodal areas, 
view factors, and thermal conductivities are given in Table IX. Steady 
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E99 2 - 34Ri 



Fig. 33. Heat fLow paths in the grid system. 30 cm P in - 200 W. 

P A = 20W. 

Acc 
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TABLE IX 


Nodal Areas, View Factors, and Thermal Conductivities 
for the JPL 20 cm Thruster 


View Factors 


F = F 
2,4 *9,11 


3, / 


10, 11 

10, 12 


12 , 12 


14, 15 
14, 17 
14, 18 


= F 15,15 
" F 15, 17 


15, 18 


16, 17 
16, 18 


18, 18 
18/17 


= 0.27 
“ 0.35 
= 0.24 
= 0. 14 
= 0.15 
= 0.29 
= 0.07 
0.06 
= 0.68 
= 0.26 
- 0.06 
= 0.50 
= 0.25 
= 1.0 


Areas, ft 2 


= 0.0351 

= 0. 0225 

- 0.0147 

= 0. 0297 

= 0.0566 

= 0.0297 

= 10 9 

= 0.105 

= 0.0675 

- 0.044 

= 0.089 

= 0.169 

= 0.089 

= 0. 175 

= 0.112 

= 0.073 

= 0. 148 

= 0. 282 
= 0.148 



Thermal Conductivities 

, Btu/hr/°F 


C 3, 1 0 

0. 299 


C 10, 16 

0.681 


C 6, 13 

2. 01 


C 13, 19 

4. 55 

r a 

C 2,3 - 

C 9, 10 = C ] 5, 36 = C 4, 5 = 

C 17 J 18 = C 5, 6 " C 12, 13 = 

C 18, 19 

= 10 6 ^ 
°C 


a Instead of calculating the very high conductances between these nodes, 
we insured that they arrive at equal temperatures by coupling them 
with very high conductances. 







state temperatures of 350°C and 322°C were measured by thermocouples 
located at the center and edge, of the screen electrode with the thruster at 
the operating point: 


Beam current = 

577 mA 

Accelerator current = 

8 mA 

Total voltage “ 

4 kV 

Mass utilization = 

69% 

Arc power ' = 

100 W 

Cathode power - 

100 w 


A set of screen and accelerator temperature distributions were 
calculated with a total power of 150 W incident on the screen for a 
Q /Q _ =1.5 distribution, for various values of interception power. 

The curves are given in Fig. 34. The case of 30 W of interception 

i i 

power (15 W into each node and N^) gives accelerator center and 
edge temperatures of 354°C and 326° C, which are seen to be in 
excellent agreement with the measured values. The value of 30 W for 
the interception power was chosen to be equal to that observed in the 
comparison thruster, the 150 W total power was assumed as a reasonable 
fraction of the total discharge power (200 W), and the value. Q /Q d =1.5 
was chosen to give best agreement between theory and experiment; this 
last value is also quite reasonable. It is felt that, within the uncertainties 
of the parametric choices, the closeness of the calculated and observed 
temperatures shows that the analysis is a useful tool. 

C. Stress Analysis of Electrodes at Operating Temperature 

In Section III-B a thermal model of the grid system was developed 
for predicting the temperatures of the electrodes under normal (ion- 
perturbed alignment) operating conditions. Since the current practice 
is to align the electrodes precisely during assembly, at room tempera- 
ture, the differences between the electrode operating temperatures and 
room temperature represent thermal-stress producing temperature 
gradients. We believe that the stresses produced as the grid system 
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undergoes heating to operating temperature produce the largest 
axisymmetric misalignments which are likely to occur, A significant 
effort has therefore been devoted to predicting these misalignments 
quantitatively. Although it is the natural next step, the extension of our 
analysis to cover misalignments resulting from nonaxi symmetric thermal 
stresses would have been too time consuming and expensive for the Phase I 
study. We will note here only that this problem is of concern for operating 
thrusters in an array. 

A computer program which analyzes shells by the finite element 
method was used to study the thermal -buckling of the electrodes. The 
computer programs can handle shells of revolution subject to axisym- 
metric or asymmetric mechanical and thermal loads. An electrode is 
divided into washer shaped regions which form the elements of the model. 
Thermal and mechanical loads are prescribed on each element and the 
computer program solves the equations of motions for the deformations. 
Since the dimensions of each element can be specified independently, 
initially dished electrodes can be considered. 

Because the stress analysis program considers solid plates, a 

literature search was conducted to determine whether an equivalent set 

of thermal and elastic parameters could be used in the solid plate 

analysis in order to simulate the actual punched plate electrode. It was 

found that such an equivalence is possible and has been used in the 

analysis. The values of equivalent elastic constants for the perforated 

5 

plates, inferred from the paper by C^Donnell and Langer, are the 
following 


Young 1 s Modulus (%/°F) Poisson Ratio 
Screen 2. 15 x 10^ 0. 7 

Accelerator 10. 8x10 CO. 36 

1. Nominal Case 

Thermal warping of the screen and accelerator electrodes 
was calculated for the temperature distributions which correspond to an 
incident screen power of 200 W, with = 1* 5, and 20 W of 
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interception power (Case 1 of Fig. 32). The radial and axial displace- 
ments of the electrodes are shown in Fig. 35. Both electrodes were 
initially dished to a depth of 0. 010 in. , and no mechanical loading is 
present. Under these conditions the small perturbations in alignment 
which occur are (1) a decrease in electrode spacing of 0. 003 in. 
maximum at the center and (2) radially outward motion of the screen 
holes with respect to the accelerator holes, the maximum being 0. 0044 in. 
at the edge. 

From this warping calculation some indication is given of the 
behavior of a grid system (initially aligned at room temperature) when 
it is heated to operating temperature. The main effects are as follows: 

a. Radial and axial displacements of both electrodes will 
occur. 

b. Both displacements will be greater for the screen 
because it becomes hotter with a larger temperature 
gradient, has a larger open area, and is thinner than 
the accelerator electrode. 

c. The misalignments which result from the above dis- 
placements are (a) a decrease in the electrode spacing, 
and (b) a transverse (radially outward) motion of the 
screen holes with respect to the accelerator holes. 

d. The magnitude of the decrease in spacing resulting 
from bowing of both electrodes is most dependent on 
the temperature gradients; the magnitudes of the radial 
diaplacements depend on the over-all temperature levels 
of the electrodes. 

2. Sensitivity to Mechanical Design 

The above warping calculations revealed the main buckling 
motions found in the grid system. It will now be shown that the magnitudes 
of these motions are highly sensitive to the grid system mechanical design. 
The importance of three design points will now be considered: 
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(a) the sensitivity of the axial displacement of an electrode (after heat- 
ing) to the depth to which it was initially dished, and (b) the influence of 
radial clamping forces such as those produced by the electrode supports, 
and (c) electrode thickness. 

a. Initial Depth — If warping of the screen electrode is cal- 

culated under the same conditions (i. e. , temperature distribution) as 
above but for an initial dish depth of 0.040 in. (instead of 0.010 in.), 
the much larger axial displacement (shown in Fig. 36) is found. To 
look more deeply into the dependence of axial displacement on initial 
depth, the axial displacement caused by heating to a given temperature 
distribution was found for electrodes of initial depths ranging from 
0. 004 in. to 4 in. Again, no mechanical loading was present. The re- 
sulting curve (for the axial displacement at the center) is shown in 
Fig. 37. The appearance of a peak in the axial displacement curve was 
not expected, and an explanation was sought. For larger values of 
initial depth the curve approaches linear thermal expansion of about 
0. 002 in. / in. , which was expected in view of the final average electrode 
temperature of 370°C and expansion coefficient of 4. 5 x /°.C. 

An explanation of the decreasing portion of the curve (Fig. 37) 
was difficult to find because of the complicated nature of the two dimen- 
sional thermal stress problem being solved and because the competition 
between buckling forces is buried in the computer program. A satis- 
factory explanation of a qualitatively similar axial displacement versus 
initial position curve (e. g. , with a peak) was found, however, for an 
analogous thermal buckling problem. This analogous problem is simple 
enough to solve analytically and is discussed in detail in Appendix II. 

The effect of the initial depth on the axial displacement was in- 
vestigated for initial depths ranging from 0.004 in. to 0. 1 67 in. This 
was done using the thermal loading of Case 1, Fig. 32, screen electrode 
with no radial constraint, and for electrode thicknesses which ranged 
from 0. 020 to 0. 100 in. Two different support configurations were con- 
sidered, both having center supports, the first having edge support at 
9 in, and the second at 12 in, diameter. The cases are summarized in 
Table X. 
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DISTANCE FROM ELECTRODE CENTER, in. 


Sensitivity of axial displacement to initial depth 






AXIAL DISPLAY 



Fig. 37. 


Axial displacement as a function of initial depth 




TABLE X 


Change in Shell Depth for Different Initial Depths 


Center Support and Edge Support at 9 

in. Diameter 



Thickness, in. 



0,. 020 

0. 035 

0. 050 

0. 100 

Initial Dish, in. 


Change 

mils 


0.004 

2.98 




0,020 

9. 15 


1.6 


0.040 

14. 3 

5. 78 

3. 01 

0. 812 

0.083 

1 5. 5 


5. 27 


0. 167 

11.5 


6. 56 


Center Support and Edge Support at 12 in. Diameter 

0.004 



0. 09 


0. 020 

5. 02 




0. 040 



0. 906 


0.083 



1. 77 


0.167 



2. 85 



Note it is generally desirable to minimize the initial dish to keep the 
axial displacement at the lowest possible value. Some dish is needed 
to insure that the resulting bowing will be in the right direction. For 
the 0. 020 in. thick case it is seen that an initial dish of « 0. 080 in. re- 
sults in a maximum axial displacement; for initial dish values above 
this the axial displacement decreases slowly. For the range of initial 
depths considered, however, we do not see this effect for the 0. 050 in. 
case. 

b. Thickness — Table X shows the effect of thickness for an 

initial dish of 0. 040 in. with a center support and edge support at 9 in. 
diameter. These data, plotted in Fig. 38, indicate the quite rapid 









increase in axial displacements for electrode thicknesses less than 
~ 0.040 in. Thus, thermomechanical considerations dictate the use of 
the maximum permissible thickness, consistent with the constraints 
imposed by discharge efficiency and ion-optical characteristics. 


c. RadiaL Clamping — The computer program used to calcu- 

late the buckling of electrodes with no mechanical loading was also 
capable of finding the buckling of a radially clamped electrode. A strik- 
ingly large dependence of axial displacement on constrained radial dis- 
placement was found. Before* we give the numerical results, we con- 
sidered a simple example of mechanical deformation of a parabolic dish 
which points out how geometrically magnified axial displacements result 
when the radial displacements of an electrode are constrained. 


Consider a dish of depth h 
a parabola of revolution: 


o 


and radius R taken in the form of 
o 


2 = h °fe) 2 - • 


Suppose that the radius of the given dish is reduced mechanically in 
such* a way that the surface area of the dish remains constant. The re- 
sult is that the dish will deform into a new shape with an increased 
depth. We wish to find the relationship between the change in radius 
AR and the change in depth Ah. This is easily done if it is assumed that 
the new shape of the dish is also a parabola of revolution with radius 
R q - AR and the original surface area. With the* original dish shallow, 
i, e. , 


h « R , 
o o 


its surface area is 


S = tt{R 2 + h 2 ) . 
o o o' 
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The equation for the surface area of the deformed dish is 

,2 , 2 . 

ir(R + h ) = S q , 

which can be solved for the new depth h as a function of the new radius 
R to obtain 

h =~\l R 2 - R 2 + h- 2 . (32) 

V o o 

In Fig. 39 the increase in depth (Ah. = h - h Q ) has been plotted as a func- 
tion of the decrease in radius (AR = R - R) for values of initial depth 

o 

equal to 0.004 in., 0.02 in., 0.04 in., and 0.06 in. The important 
effect to observe is the enormous magnification between AR and Ah. 

From (32) the derivative dh/dR is easily found to be 

dh R 

dR = “ h ' 

Thus, the geometric magnification has magnitude R/h which is indeed 
very large for a shallow dish. 

The effect of radially constraining an electrode which is thermally 
expanding was calculated by the Stress Analysis computer program. The 
case (shown in Fig. 36) of the screen electrode expanding from 0.040 in. 
initial depth to a final depth of 0. 062 in. was used. Without radial con- 
straining forces the radial expansion of the electrode is 0. 01260 in. In 
the subsequent calculations a radial boundary condition was included 
which limited the radial expansion to values between 0. 0124 and 0. 0126 
in. This radial constraint caused additional buckling of the electrode, 
as shown in Fig. 40. The curve is the final axial position at the center 
of the electrode versus the difference between free radial expansion 
(0.0126 in.) and the values specified by constraint. Thus, if clamping 
forces reduce radial expansion by 0,0001 in. the electrode center will 
move an additional 0.0135 in. (from 0. 062 to 0.0755 in. ). 
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Fig. 39. 

Geometric effect of radial clamping. 
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Consider the zeroth order approximation in which the electrode 
supports are replaced by a constraint which fixes the radial boundary at 
a given radius. The case of interest is first run with no radial con- 
straint and the outer edge expands to a certain value; the case is then 
re-run with this value minus a AR constraint value which was chosen 
to be 0. 5, 1, and 2 mils. Table XI shows the sensitivity of the axial 
deformation to this type of constraint for a variety of cases. For ex- 
ample, a radial constraint of 0, 5 mils increases by more than an order 
of magnitude the axial displacement, compared with the no constraint 
(AR = 0) case. The above data are plotted in Fig. 41. Further study 
is necessary to predict the actual equivalent radial constraint imposed 
by a given discrete support; however, it appears certain that the opti- 
mum electrode support will have the highest possible compliance in the 
radial direction. 

TABLE XI 

Change in Shell Depth for Different Radial Constraints 
(0.050 in. thick; Center and Edge Support) 



Initial Depth, in. 

0.004 

0.040 

0.083 

0. 167 

Radial Constraint, 
mils 


Change in Depth, mils 

- 

R = 0 

0. 0925 

0.906 

1.77 

2.85 

0. 5 

2.1 

15. 6 (42. 9) a 

17.9 

13.8 

1.0 

4.11 

30.4 (81.3) 

34.1, 

25. 1 

2.0 

8. 14 

60. 1 (158) 

66.7 

47.8 


rvalues in parenthesis are for 0.020 in. thick electrode 
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Fig. 41. (a) Effect of initial shape on thermal distortion, (b) Effect of radial 

damping on plate distortion. 





d. 


Temperature Distribution — The effect of temperature 
distribution is shown below for a 0.020 in, thick electrode with an initial 
dish of 0. 040 in. having a center and edge support: 

o AT 

T_, C Center to Change in 

Outer Edge, °C Shell Depth, mils 

442 35 5. 02 

293 30 4. 29 

In this case, which has no radial constraint, the deflection is 
relatively insensitive to the absolute temperature level. On the other 
hand, the magnitude of radial displacement is nearly directly propor- 
tional to the temperature gradient. However, with a radial constraint 
of a given AR, the absolute temperature level would be expected to 
exert a strong effect. 
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IV. EXPECTED SINGLE ENGINE THRUST MISALIGNMENTS // 


A. Introduction 

In this section the perturbations in the thrust vector of a single 
engine which result from a number of modes of electrode misalign- 
ment are given quantitatively. The perturbations may be changes in 
the thrust magnitude and direction, or in the location of the center of 
thrust. Each mode of misalignment is analyzed independently. The 
effects of multiple modes can then be obtained by superposition. 

B . Elementary Modes of Electrode Misalignment 

The beam intensity profile from an aligned Kaufman thruster is 
axisymmetric but varies radially from a maximum J q at the center to 
a minimum at the edge. For typical thrusters measured beam pro- 
files can be fit well by a gaussian curve of the form 


J(r) = J q e 


-(€>'•*© 


where r Q is the radius of the thruster. For our example we have 


J 

o 


= 8 


mA 

2 


cm 



mA 
2 ’ 

cm 


r =15 cm . 
o 


Since these values of J q and correspond to the densities of current 

leaving holes, the total thruster current is given by 


I 

o 



2 ir r J(r) dr 



where f^ is the opened area fraction of the accelerator. Carrying out 

the integral yields the formula for the total current I in terms of the 

_ o 

average current density J q : 
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-> J - J, 

T 2 -=• . T - o 1 

o «o o A’ J o " log J q - log 

For our thruster f A = 0.4, "j = 4.3 (mA/cm 2 ) and I. . = 1.2 A. 

The thrust vector from a thruster with misaligned electrodes 
will be found by summing the perturbed thrusts contributed from 
each hole pair. First we consider the results of the ion optical 
study which gives the current density for a single hole pair under 
three modes of misalignment: (a) change in hole spacing, (b) rela- 

tive displacement of the holes transverse to the beam, and (c) tilt 
of one hole out of its original plane (see Fig. 42). It turned out that 
only the first two modes produce significant thrust changes. Expres- 
sions for total engine thrust vector misalignment will then be derived 
for various electrode motions. 


Changing the distance d between holes causes the plasma 
sheath position to move and results in a change in the total current 
from the hole pair. To include this effect in our analysis we let the 
current density J also be a function of the hole spacing d. For aligned 
electrodes the spacing is uniform d Q (=0.090 in. ) and the current density 
is 



The dependence of the current density on spacing will be assumed 
linear and for small changes in spacing Ad 


J(r, d Q + Ad) = J(r, d Q ) jl + aAd 

where a is the percentage change in current density per unit spacing 
change. The value of a determined by the ion optical study is 6.5%/mil. 
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Fig. 42. Individual modes of aperture misalignment. 



The transverse displacement of the accelerator hole with 
respect to the screen hole by an amount s causes a deflection of the 
total thrust produced by the hole pair through an angle 0 (see Fig. 42, 
case b). The total current from the hole pair is not affected, to first 
order, and is still governed by the distribution J(r,d Q + Ad). The 
deflection angle 0 is only a function of s to first order and is linear 
for small values of s. Accordingly, with As small 0(As) = (3As 
where |3 is taken constant and equal to 0.411°/mil. The densities 
of the axial and transverse components of the currents produced by a 
hole pair with sidewise misalignment As are then given by 

J (r, d + Ad, As) = J(r, d + Ad) cos 0(As) 
z o o 

J t (r, d Q + Ad, As) = J(r, d Q + Ad) sin 0(As) 

The magnitude and direction of the thrust from a single engine 
will now be calculated for a number of modes of electrode misalign- 
ment. Superposition of the contributions from each hole pair is achieved 
by integrating the current densities and over the emitting elec- 

trode area. The misalignments considered will begin with the simplest, 
uniform spacing change, and progress to nonaxisymmetric electrode 
spacing s . 

1 . Axisymmetric Electrode Misalignments 

For axisymmetric electrode misalignments the current 

densities J (r, d + Ad, As) and J {r, d + Ad, As) are constant on 
z o t o 

rings of radius r. Different modes of misalignment can be described 
by choosing Ad and As as functions of r. 

Case 1 . Ad = constant, As = 0 
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The simplest case corresponds to uniform spacing. The total thruster 
current is axial and is given by 



I 


I 

o 


1 1 + aAd 


Ad = constant . 


The total engine thrust is T(Ad) = T q \1 + aAd) where T q , the total 
(unperturbed) engine thrust, is related to I by T q = k ?Jv 1 Q = (21.25) 
(1 . 2 A) =25.5 m-lb. 


Case 2. Ad = 0, As £ 0 

The second case corresponds to a transverse displacement of 
the accelerator electrode with respect to the screen electrode through 
a distance As. Now the hole current densities are 


J z = J(r,d Q ) cos (pAs) 

= J(r, d ) sin (pAs) 

and the total axial and transverse thrust components are 


T = T cos (pAs) 
z o r 

T^ = T q sin (pAs) 

Case 3. Ad = constant, As = rjz^ . 

The third case corresponds to rotating the accelerator grid about 
the thrust axis through an angle J>. 
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In this mode of misalignment all accelerator holes at radius 
r are rotated through an arc length x</>, where ^ is measured in 
radians. A ring of hole pairs misaligned in this mode will produce 
axial thrust and a torque about the z axis. The total axial thrust 
from the engine is 


T z (izf, Ad) = Wvf J 2-rrr J(r, d Q + Ad) cos (flr^) dr) f^ 


The total torque produced by the thruster is 


t^(/, Ad) = Wvf J 2trr^ J(r, d Q + Ad) sin (pr^) dr )• f ^ 


For the values of 4 of interest the argument f3rjz( remains small 
(less than 7°); thus the approximations 


o i 9 

cos 9 = 1 - — =- 

sin 9 =9 


maybe used to simplify the integrals. Both integrals can then be car- 
ried out and yield 


T z (^,Ad) = T q (1 + aAd) U 




T z (f(, Ad) = r Q T o (l + aAd) 



where the angle 


Pr o ^ 


is in radians . 
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4 - m ■ — 

The fourth case corresponds to a thermally warped set of electrodes 
whose shapes are parabolas of revolution. This mode of misalign- 
ment will occur in operating thrusters and has been discussed in 
Section III. The electrode spacing d(r) is the difference between 
the screen and accelerator positions z (r) and z (r). For para- 

S 3* 

bolic electrode shapes the functions z (r) and z (r) have the form 


Case 4. Ad = d 


and 


z (r) = p ( 1 - f—') jd 

s' s\ r o 

\ \ o/ / 

z (r) = d < 1 + p (1 - f-—\ ) r 

a °\ a V Vo/ J ) 

d(r) = z - z = d < 1 + {p -p) ( 1 - ) / • 

a s o vr a *s' \ lr J ) \ 


Since the screen warps more than the accelerator, p > p . We there- 

S cL 

fore take p = p - p >0 and use the above representation of d(r) in 
this case. The entire engine thrust will be axial but the magnitude will 
depend on the differential bowing which is characterized by the param- 
eter p. The total engine thrust is 


T z (p) - ^ k n/V J 27 rr J (r , d Q + Ad) dr j f ^ 

o 


where 
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After carrying out the integral we obtain the formula 

2. Nonaxisymmetric Misalignment 

It has been experimentally observed that the thermal 
deformation which occurs when electrodes are heated is not entirely 
symmetrical about the central axis. This occurs because all the sheet 
stock from which the electrodes are manufactured is formed by rolling 
and therefore has a "grain. M This exists even after thorough annealing 
and stress relieving and causes a cylindrical mode of warping for the 
electrodes. While it is not possible to predict analytically the extent 
to which this component will be present, it is desirable that the effects 
of such cylindrical curvature be documented. This will permit the pro- 
posed measurements of thrust deflection and electrode shape to be 
correlated. 

Expressions for the axial positions of the electrodes in this mode 
of misalignment are 


z = d (l + X. sin — ^ 
a o V a r Q J 

where rj is a rectangular coordinate in the direction of roll in the 
electrode material. The spacing between electrodes is then 


d = 


z 

a 


z = d 
s o 


(K 


\ \ ‘ Trn ) 

s' r ( 

o ) 
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Assuming that A. > \ , we let A. = \ -A. >0 and have 

s a s a 

?! • 

O J 

Since the electrode spacing is symmetrically distributed about d^ the 
total engine thrust is T . The thrust intensity is increased on half 
of the electrode (for ^ > 0) and decreased on the other half. The 
center of thrust is therefore shiftedtoa new off axis location (0,^). 

When written in the transverse coordinate system (£, r\ ) the 
integral for the thrust vector displacement T) o is 


T* 



where r\ = r sin v , as is shown in Fig. 43. The (i, rj } coordinate 
system is introduced for the convenience of having the thrust vector dis- 
placement along one axis (rj ) . In general, the (£,*n) system will be 
rotated by an angle from the thrust stand coordinates (x, y). If 

(Ax, Ay) denotes the location of the shifted center of thrust in (x, y) 
system, then Ax and Ay are given by 


As - sin v 
o o 


d = d \ 1 - X sin 


Ay = Ti cos v 
J 'o o 


The angular integral in the expression for T) can be carried out to 
give 


r| - - 2d r aX — ^ f u^ d * (ttu) 

° o J J 


2 o 

-u log 


du 


o o 
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where ^ Bessel function of order 1, and the substitution u r/r^ 

was made . Finally, by evaluating this last integral numerically we 
obtain 

( J 

ri = - \ ) (0. l)d r a — 

° ) ° ° J 

( o 

C. Expected Thrust Vector Perturbations 

A summary of single engine thrust vector perturbations caused 
by the five elementary modes of electrode misalignment is given in 
Table XII. Applications of these formulas to a variety of physical situa- 
tions will now be presented in order to give a feel for the implications 
of our study. 

1 „ Uniform Spacing Change 

The over-all effect of any axisymmetric spacing change 

Ad between the electrodes is a change in the (axial) thrust magnitude 

T . The relationship between T and Ad is linear and is charac- 
z ^ 

terized by the constant a = -0. 65%/mil which was determined in the ion 
optical study. Since uniform spacing changes can be easily prevented 
by the electrode supports, this mode of misalignment should not present 
3 , significant problem. For an uncertainty of 2 mils m spacing (which 
may be due to assembly tolerances) the uncertainty in thrust magnitude 
is 1.3%. 

2 . Transverse Displacement 

The transverse displacement As of the accelerator with 
respect to the screen causes an angular deflection 0 in the thrust 
vector. For small displacements the relationship between 0 and 
As is 0 = pAs where (3 is the constant 0 . 41 1 /mil which was found 
in the ion optical study. The expected amount of displacement is 
highly dependent on the particular electrode support design employed. 

The best technique for estimating As is perhaps by vibration testing. 

In any case, uncertainties of 10 mils (for example) relative in trans- 
verse positions of the electrodes cause uncertainties of 4° in thrust angle. 











3. 


Rotation About Thrust Axis 


A rotation ^ of the accelerator with respect to the 

screen about the thrust axis causes a decrease in axial thrust and 

a torque t about the thrust axis. The expected amounts of rotation 

also depend highly on the electrode support design. Upon substitution 

into the formula for t we find that 

z 


t ,.(4.3**10- s i£r) ^ • 

I 

Thus, 1° of relative electrode rotation will cause a roll axis torque of 
4. 32 x 10^ in. - lb . 

4. Parabolic Electrode Spacing 

The change in thrust magnitude has been calculated for 
the parabolic electrode spacing shown below 


E 1010-27 



This mode of misalignment results from the thermal buckling of the 

electrodes. When the electrode gap decreases from d Q at the edge to 

n • d at the center (p < 1), the change in thrust is 
^ o 


a P d o 



0.617 apd 
r o 



pd 


o 


Thus if p = l/2, pd Q = 45 mils and the change in thrust is 18%. 

As was discussed in Section III-B the expected amount of buck- 
ling, which we have characterized by the parameter p, is highly 
dependent on the initial shape of the electrodes and the amount of radial 
clamping by the supports. Since the dependence of p on initial shape 
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is nonlinear as shown in Fig. 37, and the sensitivity of p to radial con- 
straint is severe (see Fig. 39), the uncertainty in knowing p and 
therefore the thrust magnitude is large. 

5 . Cylindrical Buckling 

A particular cylindrical non- axisymmetric mode of 
electrode misalignment was considered because it has been observed in 
operating thrusters and because it apparently cannot be eliminated. The 
"grain” induced in the electrode material during rolling is the cause of 
this buckling mode. The electrode spacing is shown below 



dOl) = d Q 1 


X. sin 


TT~n 

r 

o 


as a function of tj , the direction of roll. In this mode of misalignment 
the thrust is axial and has magnitude T^, but the center of thrust is 
displaced by an amount t| q . The relationship between X and r\ Q is 

tt = (0 . 65 in. ) X . 

Thus, for Xd Q = 20 mils, X = (2/9) and t\ q = 0.145 in. 
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V. 


THRUST STAND 



In order to carry out the experimental verification phase of this 

program, a device is required to measure the location of the thrust 

vector from the thruster described above to within ± l/4° of angle and 

to determine whether any couple exists about the thrust axis to an 

-4 

accuracy of 5 x 10 in. -lb. Competitive methods of making this 
measurement were analyzed. Conventional techniques for measuring 
either thrust or ion trajectories are investigated first. Each of these 
(viz. , scanning Faraday cups, torsion wire suspension, and ion beam 
collection reaction) was eliminated for this application. At this point 
two new concepts of thrust measuring system were considered. It 
appears that each has sufficient sensitivity to measure the thrust vector 
directly, but that one may possess some inherent advantages. The 
systems considered are discussed below. 


A. Scanning Faraday Cup 

The basic concept here is to use a Faraday cup to establish 
the number and direction of the ions at each point in the beam 
and from this information to infer the thrust direction. As with all the 
methods compared the limiting accuracy must be ± l/4°. A fundamental 
problem associated with the various ways in which this general technique 
may be applied is that the ions in the beam do not follow laminar tra- 
jectories. For this reason, it is necessary actually to measure the 
direction that the ions are traveling at each point. 

The most direct way to accomplish the above is to make the mea- 
surements at a point sufficiently far downstream that the actual source 


The substantial contributions of Mr. S. Hansen of HRL to the concept 
and analysis of these thrust stands are gratefully acknowledged. 
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may be approximated by a point source. With the above required 
accuracy, the source must subtend an angle less than l/E° at the cup. 
For a 30 cm diameter thruster this gives a beam path length well over 
100 ft. Not only are facilities of this size not available, but scattering 
in the beam would make the results uncertain. 

The second possibility is to employ a Faraday cup relatively 

close to the thruster. The cup must now be designed to accept only 

ions which have emanated from an area on the source which subtends 

a 1/2° angle or less at the collector. In other words, the cup must be 

preceded by a very narrow angle collimator. In principle at least, this 

device may now be used to measure experimentally the beam intensity 

as a function of angle (between Faraday cup beam axis) at each point in 

the beam. This measurement technique of course would generate an 

enormous amount of data. The experiment is difficult not only because 

of the accuracy required to position and control the location of the cup, 

- 8 

but also because of the low signal level available (~ 10 A) after the 
selective collimation. Actually reconstructing the thrust vector from 
such measurements would no doubt require considerable computer 
analysis . 

Another important factor which mitigates against this method is 
that collection of the data to establish a single thrust vector might re- 
quire several hours, a period which may be much longer than the 
stability of either the measuring equipment or the thruster. 

In summary, it does not appear practical to construct the thrust 
vector to the desired accuracy from detailed measurements of ion 
trajectories in the beam. The principal difficulties arise from the very 
large amount of information required to define the beam with sufficient 
accuracy. 

B. Torsion Wire Systems 

The basic concept here is to hang the thruster from a thin 
wire and then measure the angle through which the wire twists 
because of the torque created by unbalanced forces. A number of 
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variations of this general idea, such as suspending the thruster between 
two wires (as in a taut wire galvinometer) and using nulling techniques 
to maintain the net angular deflection at zero, all fall into this same 
general category. 

With this type of suspension system it is conceptually possible 

-4 

to measure torques of the magnitude of interest here (5 x 10 in. -lb). 
However, there are a number of practical difficulties which seriously 
limit the application of the technique. 

• Electrical Leads — A number of electrical connections 
must be made to the thruster. Each of these will stiffen the suspension 
and reduce sensitivity. A possible solution is to allow each lead to dip 
into a separate mercury pool, thus providing a friction free contact. 

It is very difficult to estimate the chance of success of this technique; 
it is at best a very cumbersome one. 

• Stability — This arrangement is very susceptible to 
pendulum type instabilities. In principle, these can be damped out, 
but only with a sacrifice in over-all system response. 

• Nulling Mechanism — Considering the device as a 
simple pendulum, it may be shown that the engine thrust (25 mlb) will 
cause motions of the thruster system of the same order of magnitude 
as those created by the torsional forces about the wire axis which we 
are trying to measure. This system thus requires a relatively sophisti- 
cated nulling system which automatically balances the main thrust 
vector at all times while at the same time measuring the much smaller 
(by a factor of approximately 100) force necessary to compensate for 
the thrust variation which we wish to measure. 

• Interpretation of Data — A system such as this would be 
designed to measure only the torque about the axis of the suspension 
wire. Inherent in the system are a number of other degrees of freedom 
which may be excited by perturbations in the thrust direction. Coupling 
between these may make interpretation of the experimental results very 
difficult. 
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In summary, while this type of system conceptually has just suf- 

_4 

ficient sensitivity to measure the torque specified (5 x 10 in. -lb), a 
number of practical problems make application of the basic principle 
difficult, if not impossible. 

C. Collectors 

Under ideal circumstances it is possible to measure 
thrust produced by an ion thruster by measuring the effect of the 
impinging ions on a collector. The principal practical difficulties 
which arise are associated with the dissipation of the beam power 
(~ 2 kW must be carried away from the collector by some type of 
forced cooling system, usually water). It is impractical to consider 
water leads with sufficient flexibility to permit the collector to respond 
to thrust levels of the magnitude of interest here. Therefore, it is 
not possible to measure the thrust direction by observing the reaction 
of the collector to the thrust. 

D. Stiff Column Thrust Stand 

The heart of this design is an extremely sensitive tilt- 

meter recently developed at HRL which is capable of resolving 

_q 

deviations of 10 rad from the horizontal. The basic concept is 
illustrated in Fig. 44. The thruster is so positioned that both the 
gravitational vector and the nominal thrust vector are parallel and 
lie along the axis of a slender hollow column which both supports the 
test package and contains the necessary electrical leads to operate 
the thruster. The base of the column is supported on a massive 
platform stabilized by a servo loop which is controlled by a tiltmeter. 

A second tiltmeter mounted on the platform supporting the thruster 
monitors small angular fluctuations which occur when variations in 
thrust direction cause the slender column to deflect. 

The sensitivity of this system may be estimated as follows. 

The angular deflection of a slender column is given by 
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rad 


6 


1 Ft 2 
' I El 


where 

F = force perpendicular to column 

£ = length of column 

E = modules of elasticity 

•I = moment of inertia of cross section of column. 

For a 6 in. long quartz rod, l/4 in. in diameter, used in the proposal, 
-2 -8 

0 = 10 F. With the sensitivity of 10 rad specified for the tiltmeter, 
_ 5 

forces of 10 lb may be detected. 

The basic concept of monitoring the position of a flat "pancake" 
shaped bubble formed in a liquid under a flat plate, which is the heart 
of the tiltmeter design, may be employed in two ways to achieve the 
desired result here. The first (Fig. 45) is to provide a servo system 
to maintain the platform on which the thruster is mounted in a level 
position using the bubble as an accurate null indicator. If the servo 
loop has been precalibrated it may be used to directly read out the 
force necessary to counteract that resulting from thrust deviations. 

The second (Fig. 46) technique is to use a complete tiltmeter which 
includes not only the bubble sensor but the necessary suspension and 
servo system to maintain the bubble in a level position on a moving 
platform and to indicate table declination directly. 

A detailed design of a representative system of this type was 
presented in the first quarterly report under Phase I of this contract. 

It was demonstrated that a 12 in. long stainless steel tube, 3/4 in. o. d. 
with a 0. 020 in. wall would serve satisfactorily to support a 32 lb load 
while providing the required flexibility for use with a standard HAC 
4 in. diameter tiltmeter. A mathematical analog was developed in the 
above quarterly report which indicated that with typical values a 19 mV 
signal would be available at the minimum thrust deviation of l/4°. The 
frequency response can be made essentially linear up to 1 cps. 
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Figure 47 depicts the perturbing forces associated with the 
misalignment of the thrust vector and the center of mass. Ideally, 
these forces should act along the axis of the supporting column. Dis- 
placements of these two forces from the ideal position are represented 
by x and y, respectively, in Fig. 47. The column of length t is 
supported on an accurately level table and makes an angle (3 with the 
horizontal. Assume first that |3 = 90° and that the two perturbing 
torques (taken about point A at the base of the column) are 

Lj = (T cos a)(x) 

L 2 = Mg y 

where 

M = mass of the system (assume 25 lb) 

T = total thrust (assume 0. 025 lb) 
a = thrust deviation (> l/4°) . 

The torque produced by the thrust component that we wish to measure 
is 


L 3 = T(sin a) l 

where & = column length (assume 12 in.). To assure that 
and we must have 


x < 0. 050 in. 


and 


y 


<5x10 


-5 . 

m. 
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These calculations show that there is a relatively large tolerance on 
the alignment of the thrust vector with the axis of the supporting 
column, while very little shift in center of mass position is permis- 
sible. 

Consider first the thrust vector misalignment. For a perfectly 
symmetrical thruster the thrust vector should pass through the exact 
center of the ion optical system. This can certainly be specified to 
better than the 0. 050 in. permissible error calculated above. Deter- 
mining the probability that a shift of this magnitude will occur during 
testing should be one of the goals of the analytical study. It is fortunate 
that reasonable lateral shifts in the thrust vector location are tolerable, 
because it is not possible to separate the two torques 1^ and by 

varying the total thrust, since each varies directly with the thrust 
magnitude and angle. The only other variable which even in principle 
permits the separation of the two effects is the misalignment x. As 
shown below, this cannot be adjusted during a run because of the very 
large effects associated with shifting the position of the center of mass. 

The force on the thrust table due to gravity acting on the center 
of mass of the table and its load (the thruster, etc.) must pass accurately 
through the pivot point A at the base of the support column. Since 
neither the center of mass nor the pivot point can be located to the 
necessary accuracy, a technique for mechanically balancing the system 
is required. The location of the thrust vector physically defines where 
the thruster must be located with respect to point A; therefore, balancing 
must be done by the addition of weights to the table, possibly using one 
of the many techniques developed for modern analytical balances. 

Final precision balancing may be done with automatic nulling incorpo- 
rated into the system to maintain the bubble level during the test. 

Small shifts in the center of mass may occur during operation 
because material is being removed from or deposited on the ion 
source or its supporting structure. Other center of mass shifts may 
inevitably occur when the electrodes are moved with respect to each 
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other as part of the experimental program* These latter effects can 
be calculated and also measured directly by shifting the electrodes 
with no beam extraction and thus may be subtracted from the observed 
thrust deviation to give the true reading. Because the electrode mass 
and the amount of motion will be known, it is also possible to compen- 
sate for the center of mass motion by moving an equivalent weight in 
the opposite direction. The permissible tolerance for angle (3 may 
be directly calculated once the column length i is fixed. 

E. Floating Suspension Thrust Stand 

This new concept of thrust stand uses a suspension 

mechanism which permits the three components of the thrust and the 

❖ 

torques about the three mutually perpendicular axes to be measured 
at the same time. The thruster platform is floated in a mercury pool 
and is free to drift horizontally, or to lean away from the vertical, 
depending on forces applied. It is also possible for the system to 
rotate about its vertical axis if a couple exists about this axis. As 
discussed below, the suspension.must be carefully designed to assure 
stability and to assure that the suspension itself introduces minimum 
perturbations into the system. 

The operation of the system may be understood by referring to 
Fig. 48,. which shows a two dimensional sketch of the essential elements 
of the device, and Fig. 49, which illustrates the two types of thrust 
misalignment which may occur. The mercury surface remains 
horizontal, thus removing the necessity of placing the entire apparatus 
on a stable table, as in the previous design. In operation the table is 
releveled by the force motors F , using the tiltmeter as a reference, 
while horizontal and rotational motion is arrested by the force motors 
F and F^, respectively. The position monitors which establish these 
null positions are not shown. 


The total thrust is calculated from electrical measurements. 
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CASE 7 CASE 2 

Fig. 49. Two types of thrust misalignment. 


Consider first case 1 of Fig. 49- Here the thrust vector is at 
an angle 0 with the desired direction along the centerline. It is appar- 
ent that such a force tends both to tilt the platform and to move it hori- 
zontally. To resolve such a force, force motors F level the platform 
and the force motor F compensates for the horizontal drift. Next, 

X 

consider case 2 of the same figure, which illustrates a situation in 
which the thrust vector has moved laterally while remaining parallel 
to the centerline. Notice that if l ^ = l the same torque is produced 
about the center of rotation, as in case 1. However, there is no lateral 
force in the second case. 

The primary advantages of the floating suspension are (1) that 
it resolves the two possible types of thrust misalignment illustrated 
in Fig. 49 and (2) that a single experimental setup measures both 
angular deflections of the thrust vector and torques set up about the 
centerline or nominal thrust axis. Figure 50 illustrates the proposed 
suspension mechanism which provides the additional degrees of 
freedom. 

It is desirable to design the floating suspension system so that 
the buoyant force of the mercury supports the platform but does not 
tend to force it to maintain a particular position. In other words, by 
proper mechanical design the system may be made critically stable so 
that the buoyant forces tend neither to right it or to upset it. The 
design criteria which define such a neutral system were described in 
an earlier quarterly report. 

The preliminary design of the system incorporated a bubble 
sensor, such as is used as the sensing mechanism on the HAC 
tiltmeter, to sense the table position with respect to the horizontal. 

It was subsequently recognized that it is not necessary to maintain 
the table accurately horizontal at the null position when the neutrally 
stable suspension described above is achieved. By using linear 
position sensors referenced to the platform supporting the mercury 
reservoir to monitor the motion of the floating table in the vertical (z) 
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as well as the horizontal (x, y) plane, the dynamic control system is 
greatly simplified and made unconditionally stable. Both systems are 
described and compared in the following discussion. 

As shown in Fig. 51, the several components of thrust and 
torque are balanced by a combination of eight force motors; four are 
aligned with the vertical axis z, positioned in a plane symmetrically, 
and four are aligned in the horizontal plane. It will be demonstrated 
below that this arrangement is adequate to determine the thrust de- 
flection and torque by algebraic manipulation of the static force motor 
readings. 

The vertical force motors are intended only to balance torque 
about an axis in the horizontal plane (x, y). The main thrust (along 
the z-direction — vertical) is obtained from measurement of the ion 
engine beam current and voltage; thus z-motion is unimportant. 
Therefore, these four force motors are operated as series pairs such 
that the two always produce an equal and opposite force. This gives 
the required torque with no net z force. 

The horizontal force motors must balance torque about the 
vertical (z) axis and translational thrust components in the x-y plane. 
Therefore, the opposite horizontal force motors cannot be paired to 
always produce an equal and opposite force because they must produce 
a net balancing force in addition to the torque. 

Figure 52 shows the transducer and force motor control 
arrangement for the five degrees of freedom of interest. The (3 

x 

an d Py rotations are sensed by the (two-axis) bubble sensor or by 
four linear transducers. These signals drive the four vertical 
motors in matched pairs. The rotation about the vertical axis and 
the translation in the x-y plane cannot be sensed independently. This 
coupled motion is observed by four photoelectric transducers, which 
sense the net x-y motion of a point on the edge of the platform. Hence, 
the displacement which each horizontal transducer experiences is the 

sum of x-y center of gravity translation, plus rotation (a sin p ). 

2 
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Fig. 51. Thrust stand schematic layout 






















Each of the horizontal force motors is controlled by its associated trans- 
ducer, independently of the other three. The total system thus has six 
control loops. 

1 . Static Analysis '" 

The principle of thrust stand operation is that under steady 
static conditions, the static forces applied by the force motors exactly can- 
cel out the applied constant force of the thrust vector plus an associated 
torque. From a knowledge of the force motor readings, the thrust vector 
deviation from the normal plus the torque may be determined. 

The thrust vector deflection is determined by balancing the 
horizontal force motors. The thrust vector is represented by 


F = Fi+Fj+Fk 
x y z 

where i, j, k are unit vectors in the x, y, z directions. The angu- 
lar deflection 0 which we seek may be expressed as 


sin 0 



+ F 


2 

y 


(33) 


In order that the x and y components will balance exactly, 


F 


iy 



F 

y 


and 


2x 


+ F 


4x 


= F 


The amplitude of the total thrust F is measured independently from 
beam current and voltage. Substituting into (33), 


^The analysis of this thrust stand was carried out by J. D. Smith, 
D. Berman, and S. Harrison. 
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“*» - T [( F 2* + O' + ( F ly + F 3y) 2 ] 


(34) 


Therefore, we obtain the deflection 0 from calculation involv- 
ing only the basic horizontal force motors. This is true in theory; in 
practice, success will depend on the accuracy of the readings, and the 
required deflection resolution. For l/4°, 

sin 0 ~ Q = y ( y gg ) ~ 0* 004 rad. 

Assuming no torque about z and equal x-y components (F^ x = F^ x = 
Fj = F ^y =’ F^), £h e force measuring accuracy required to measure 
this deflection, for a nominal thrust of 0. 025 lb, is 


eF = '\( 2F -y f + ( 2F xy) 2 = F : 


Evaluating, 


F - = 0.35 xlCf 4 = 35 Ub. 

2 yg 

The applied engine torque is obtained by balancing the total 
applied and induced torques against the combined torsional reaction 
of the force motor array* The induced torque is F x R, where 
R = Axi + Ayj + Azk = the vector distance between the point of thrust 
application, and the total system center of mass. The induced 
torque components about the x, y, and z axis are, respectively, 
(AzF y - AyF z ), (AxF z - AzF x ), and (AyF^ - AxF y ). 

The torque balances about each of the three axes are 


(F 1 Z - F 3z>*H' - (F 2x + F 4x>‘v = AzF y - A V F y • < 35 > 

t F 2* - F 4 ,)'h' - < F ly + F 3y> l v = ^ - A * F * < 36 > 
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(37) 


~(F 


ly 



H 


< F 2x - F 4x>*H 


AyF - AxF - T 
J x y z 


where T is the engine applied torque, and l is the lever arm of 
z 

the force motors about the system center of gravity (Fig. 53). 

Equation (37), balance of torques about the z-axis, yields the 
applied torque T , provided we know Ax and Ay, the lateral dis- 
placement of the applied force. These are available from (35) and (36). 
The distance Az is not unknown, since it represents the distance 
between the accelerator electrode of the thruster and the system center 
of gravity. Hence 


T 

z 



0 


+ { AxF 



-4 

The accuracy required for torque measurement is 5x10 in. -lb. 
Neglecting the thrust vector coupling, the last term disappears, and 

F ly * F 2*’ etc ' Thus ' 


T = 4F i_. 
z xy H 


From Fig. 53, the lever arm length ^ « 8 in. 


xy 


5x10 

4x8 


-4 


= 1. 5 x 10 


-5 


Evaluating, 
15 p.lb. 


Thus, the applied torque measurement accuracy requirement is 
slightly more severe than that of the thrust vector deflection 
measurement. 

State-of-the-art thrust measurements down to 10 plb have 
been demonstrated in a one-axis thrust platform developed by HRL 

vU 

for the U.S. Air Force. However the six degree of freedom plat- 
form presents new problems with transducer characteristics, 
mechanical balancing, and dynamic stability. The dynamics are 
discussed below. 

Contract AF 04(61 1) - 10536. 
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2 . 


Dynamic Model 


In order to study the thrust platform dynamics, the motion 
was limited to two degrees of freedom: one translational (x) and one ro- 
tational ((3). The rotational motion is about an axis perpendicular to the 
translation axis, making the problem essentially two-dimensional* 

These two components of motion are sufficient to demonstrate 
the different transducer loop dynamics and the coupling effects between 
the various degrees of motion. Thus, this limited model can provide 
adequate insight into the dynamic characteristics of the platform, al- 
though it should be kept in mind that the real problem has thrice the 
complexity represented here. 

Figure 54 shows the schematic for the dynamic model. 

Three force motors balance the applied force components F , 

F^. The two force motors F^, F^ provide a restoring torque to 
balance the applied torque AxF + AzF + (f^F^), where Ax, Az are 
the applied force coordinates with respect to the system center of 
mass. The force motor F^ balances the horizontal applied force 
component' which we wish to measure. 

Because we wish to study the platform dynamics, we must now 
consider other forces which were not important in the static case. 

These are displacement and damping forces which exist when the sys- 
tem is not in a null condition. These additional forces consist of 
viscous damping in the mercury liquid, spring displacement and 
damping forces from the electrical leads, and buoyancy torques due 
to displacement of the float metacenter. Thus, each degree of motion 
is described by a damped spring -mass system, although some of the 
spring and damping forces may be very small. 

In block diagram notation 


S forces 


J_ 

m 

2 

s(s + S ) + 03 
x x' X 


displacement (x) 
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nsional dynamic model 






















where 

m = system mass 

I = system moment of inertia about horizontal axis 
^x ^ natural frequency, x-motion 
T x = l/s s time constant, x-motion 

Up = natural frequency, (3-motion 
t |3 = - time constant, p -motion. 

To complete the model, we also include the coupling forces 
which can cause one degree of motion to affect another degree. This 
coupling is of two types: 

1. physical coupling 

2. induced coupling. 

The first of these can be visualized by considering the translational 
displacement of a body attached to a spring not at the center of mass. 
The displacement thus causes a reactive torque, in addition to the 
translational restoring force. The second type of coupling is that 
induced by the force motors. For example, when the x-translational 
force motor acts it will induce a torque. 

The coupling transfer frunction may be written 



Force 

or 

Torque 
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3 . 


Control System Analysis and Design 


The engine thrust measuring instrument is a precision balance 
of six components of force and torque . Before engine thrust is measured the 
system must be balanced by adjusting the engine and counterweights to re- 
duce the effects of environmental vibrations which might result from vacuum 
pumps, seismic disturbances, etc. When the balance adjustments have been 
made and the mechanical zero established, the mechanical system is then 
calibrated with the electrical instrumentation and control system. 

The following physical equations describe a single degree of 
translational dynamics of the six-degree thrust platform during an 
engine thrusting condition. These equations are based on the assump- 
tion that the platform has been previously balanced in an engine non- 
thrust condition such that platform disturbance from background vibra- 
tion is completely balanced out. 

(Fg - F.^) ~ M'x + Bx + Kx 

= engine thrust force (lb) or torque 

= platform restoring force (lb) or torque 

w 2 

= total mass of system {lb-sec /in.) or moment 
of inertia 

= velocity damping (lb-sec/in. ) 

— spring rate (lb/in. ) 

= horizontal acceleration 
= horizontal velocity 
= horizontal displacement 
Linearizing and performing a LaPlace transform on the force balance 
equation yields the linear operator expression 

Ax (S) _ 1 

AF(S) (MS^ + BS + K) 


where 



M 


B 

K 

x 

X 

X 
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The thrust platform is designed to have a minimum mechanical spring 
rate K in the direction of thrust measurement. However, all residual 
spring forces cannot be completely eliminated. As the value of mechan- 
ical spring rate K approaches zero, the transfer function of the mechan- 
ical device approaches the' following 

Ax(S) s 1/B 
AF(S) S( Tl S + 1) 

where 

Tj = M/B (sec). 

It is clearly shown here' that as the mechanical spring rate (resulting 
from wire, cable', buoyancy effect, etc. ) approaches zero, the device 
assumes the characteristics of a mechanical integrator with a cascaded 
first order lag with a time constant (sec). 

In order. to achieve moderate speed operation the platform 
mechanism control loop is closed to meet the following requirements: 

1. The goal of the dynamic response of the control system 
with a 20 lb engine weight load on the platform will be 
approximately 1 cps at all force ranges. 

2. The goal threshold of the thrust stand is less than 

3. 0 (Jilb of thrust. 

3. The goal steady state error of the thrust stand is less 
than 3% of the upper limit of the 100 filb range setting. 

4. The system is optimally damped for all ranges of 
thrust and weight loads and stable with a 20 lb engine 
weight load at all force ranges. 

Since the thrust platform is basically an integrator device, the 
platform will continue to displace as long as a force imbalance exists 
across the platform within the mechanical limits of the mechanism. 

To provide a restoring force feedback to servo set the platform 
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to its steady state null position during engine thrusting, it is necessary 

to provide a voltage integrator to generate a restoring voltage E^ to 

the restoring motor that is a function of the time integral of the dis- . 

placement away from the null position. Displacement of the platform 

Ax causes an error voltage AE x to be generated in the transducer 

output. This error voltage is integrated to contribute to the restoring 

voltage E to the force motor. When the horizontal forces across 
-K. 

the platform are balanced, the platform is restored to the null position. 
The addition of a voltage or mechanical integrator (i. e. , a bubble 
sensor) into the control loop results in additional control loop phase 
shift and system instability. Therefore, in order to provide a response 
with a high degree of stability for this type of system, it becomes 
necessary to introduce a proportional plus derivative compensating 
network in parallel with the voltage integrator and in tandem with the 
mechanical integrator. 

Horizontal loop — Consider now the transfer function block 
diagram of the horizontal force balance control loop (Fig. 55). The 
open loop transfer function for the linear system is 

K (S 2 + a S + a ) 

<^iTT/n\ U 1 Lj 

s 2 + (B/M)S + 
rad/sec 




cij = Kp/Kp (rad/sec) 
a 2 = Kj/K d (rad/sec) 2 
P = K A K M (lb/V) 


Kp = proportional gain (volts/volts) 
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Fig. 55. Horizontal force balance control loop system block diagram. 




= derivative gain (volt -sec /volt) 

Kj = integral gain (volt/sec -volt) 

P = force motor power gain (lb/volt) 

6E 

— g — = displacement transducer gain (V/in. ) . 
x 

If the network gains can be assumed to be adjustable, then can be 
set approximately to equal {K/M) and a. can be adjusted to approach 
(B/M) experimentally (e. g. , the network gains^ can be set to the follow- 
ing analytically determined values). 

a x = K p /K D = B/M * 5. 0 (rad/sec) 

a 2 = Kj/Kj-j = K^/M » 1. 0 (rad/sec) 2 

If K d = 1,0 sec/rad, then Kj = 1.0 rad/sec and Kp = 5. 0. The roots 
of the polynomial are thus 

K 

s 2 + (B/M)S + K/M = s 2 + (K^/K^JS + — 1- = o 

XJ 


cr = 




- 4 


= - 2. 5 ±~ 4 


o' 1 = - 0. 2, cr^ = - 4. 8 

Oi 

Gains less than unity (1) imply passive networks." 
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The open loop transfer function is thus 

F r (S) K o (S + 0-jHS + «r 2 ) 

F s S(S 2 + (B/M)S + K/M) 

where F^,'= F^, - F R = error signal, shown in Fig. 56. Because the 
roots of the network are adjusted to equal the roots of the mechanical 
system, the open loop transfer function is reduced to 



The closed loop transfer function now becomes 



K 

o 

S + K 

o 


1 

tS + 1 


The over-all loop response acts as a first order lag with a time con- 
stant (t) equal to 1/K q , i. e. , 


Here 


•*■(*)- 
° K 


35J 


(sec) 


M = 0. 065 lb sec 2 /in 


Kj^ = 1.0 sec/rad 


6E 


x 


X 

p 


10 4 V/in. 


io " 5 ib/y. 
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Therefore, 


t ~ 0. 65 sec. 

Improved frequency response can be obtained by increasing the propor- 
tional network gain Kp, For example, let Kp = 50; the network roots 
are then 

0- = ± i V2500 - 4 

or = - 25 ± <£ V2496 
cr^ ~ 0, 02; ~ 50 

The open loop transfer function now becomes 

F r K q (S + 0, 02) (S + 50) 

FT = S (S + 4. 8) (S + 0. T) ' 

2j 

A plot of the root locus of the open loop transfer function is shown in 
Fig. 56. The percent overshoot when the system is subjected to a • 
step input can be estimated easily by referring to Fig. 57. For a 
damping factor £ = 0. 7 the overshoot is near zero. For higher values 
of damping factor it is noted that the system is heavily damped. For 
lower values of damping, the system is underdamped. For operation 
at a damping factor of 0. 7, the over-all system natural frequency 
(co n ) is approximately 65 rad/sec or 10 Hz, as shown in Fig. 57. 

This estimate will be verified in subsequent analog solutions. 

co =65 rad/sec ~ 10 Hz 

n ' 

| = 0. 7 
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The closed loop function of the horizontal force balance system is 


where H = 1. 



G 

1 + GH 


(S) 


The frequency response characteristic shows that the theoretical 
system bandwidth is approximately 20 Hz for a system gain of 
K q = 85 rad/sec. The system equations indicate that the bandwidth of 
the horizontal force system can be extended by increasing the system 
gain K. Thus, the power gain of the control loop and the mass of the 
system will determine the. system bandwidth; for a 10 Hz system 
natural frequency with a 20 lb load on the thrust stand, the open loop 
gain of the system should be at least 


K = 
o 



85 rad/ sec. 


Rotational Deflection Control Loop — Consider now the transfer 
function block diagram of the rotational force balance control loop 
shown in Fig. 58. The basic components of the system are the 
mechanical system, bubble sensor, compensating network, and power 
amplifier and electromagnetic force motor. 

In theory, any static moment resulting from an off-axis engine 
thrust component will be counterbalanced by an equal and opposite 
torque from the force motor, provided that the control loop is properly 
designed. The design of a compensating network for the rotational 
control loop is somewhat more critical, however, if a bubble level 
sensor is used.because of the dynamic characteristics of this com- 
ponent. 

Consider the basic system problem. The Laplace transfer 
function of the mechanical system in the rotational axis, neglecting 
for the moment the cross coupling effects, may be described as 
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Fig. 58. Rotational force balance control loop system block diagram. 



A0 


(S) 


1/1 


(38) 


AF, 




2£wpS + o>p) 


where 




- f(I, 


K efP 


As the effective spring rate approaches a |small value, as in the 

case of a neutrally balanced system with negligible spring rate coupling 
of the electrical cabling, the platform system dynamics approach the 
characteristics of a mechanical integrator for small force perturbations; 
e. g. , as ojp -*■ small value, 


af e “ STS"+"2'|w) * (39) 

The bubble sensor also has the characteristics of a mechanical inte- 
grator, i. e. , 


AE _ K TB 

at {S) “ s(s + s TB ) * 


(40) 


The cascading of the two devices now presents a relatively severe 
stability problem, as shown below. 


AE _ K TB^ 

AF “ 2 “ 

E S <i (S + 2ea, p )(S + S TB ) 


(41) 


The double integration term (S ) results in 180° phase shift of the input 
signal. The compensation network therefore must correct the total 
system phase shift to less than 180° if over-all loop stability is to be 
achieved. 


The first step in considering the compensating network is to 
consider a single first order lead network coupled to the output of the 
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bubble sensor output. This network would be comprised of a propor- 
tional plus derivative network: 

= k p + k d s = K P < T C S + 11 • (42> 

The effects of compensating for this phase shift with single or multiple 
lead networks were considered in detail both analytically and by analog 
simulation in the second quarterly report of Phasel of this contract, to 
which the reader is referred for details. The final conclusion drawn is 
that it is possible to stabilize such a system. It was also concluded, 
however, that a much simpler and equally sensitive system could be 
constructed by using linear displacement sensors placed at the periphery 
of the floating mass to achieve the same result. The final design is 
based on this latter concept; the dynamic analysis is then similar to 
that for the horizontal loop. 

The cross coupling effects between loops were briefly studied 
on the analog computer. The result was that if each loop was inde- 
pendently stabilized the coupling introduces no new instability. 

The basic conclusions to be drawn from the study are that 

a. the bubble sensor may be replaced with linear trans- 
ducers and that an over-all simplification of the 
system will result 

b. ' the linear transducers should be designed to provide 

maximum possible sensitivity (10^ V/in. , if pos- 
sible) 

c. . the floating mass should be as light as possible; this 

suggests exchanging the position of the forcemotor 
magnets (now mounted on the float) and the coils if 
suitably flexible leads can be demonstrated. 
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VI. THRUSTER DESIGN 

A. Introduction 

The thruster -is representative of the state of the art and designed 
so that the electrodes can be moved relative to' each other in a controlled 
manner in order to simulate misalignments which may occur naturally 
during operation. The basic requirements of. the thrust stand are that 
it be capable of resolving the thrust direction to ±1/4° of angle and that 
torques about the thrust axis of 5 x 10“4 in. -lb or greater be measur- 
able. The detailed design of the various components and the way in 
which they will be incorporated into the total experimental system are 
discussed below. 

B. Thruster Design 

It is particularly important that the thruster designed under this 
task be representative of the current state of the art in terms of per- 
formance and fabrication techniques. With this in mind, the following 
design criteria were established. 

a. 30 cm diameter beam - based on current performance, 

a thruster of this size represents a power level of approx- 
imately 2 kW. Recent system studies indicate that this 
is a useful module size. 

b. Thermionic cathode - this cathode type was chosen to be 
consistent with experimental studies at JPE. 

c. Lightweight construction - sheet-metal techniques are 
used to hold the thruster weight to approximately 11 lb. 

In order to minimize the weight, the use of aluminum (versus 
stainless steel) was considered for the discharge chamber shell, anode, 
and other structural members. Stainless steel was chosen because 
(a) it is known to resist attack by mercury; (b) it is readily welded to 
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the soft iron end plates, while aluminum must be rivetted; (c) no weight 
penalty is associated with the use of stainless steel because it may be 
used in very thin sections. 

The estimated weights of the various* components are listed in 
Table XIII. 

d. Ion optical array - a "conventional" ion optical array 

consisting of a matrix of round apertures in a' hexagonal 
» 

close packed pattern was chosen. The following dimen- 
sions were used: 

screen: 70% open area 

0. 187 in. holes on 0. 210 in. centers 
0. 030 in. thick 
molybdenum or alloy 

accelerator: apertures 0. 144 in. (based on in- 

house HRL and JPL studies) 

maximum thickness 0. 110 in. 
based on required lO^hour lifetime 

molybdenum or alloy 

spacing: nominally 0. 090 in. , but can be 

reduced to increase perveance 

a center support is proposed to 
improve stability 

Basic design - the basic design is a scaled version of the 
SERT-II thruster which was developed recently at NASA- 
LeRC. The most significant deviation from direct scaling 
is that the discharge chamber length has .been held ap- 
proximately constant, as suggested by a HRL scaling 
study. 
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Weight, lb 


Part 

Screen Collar 
Lower Ring 
Outer Shell 
Anode 

Screen Pole 
Piece 

Screen 

Electrode 

Accel Electrode 

Back Plate 

Cathode Pole 
Piece 

Magnets 

Insulator 

Assemblies 

Subtotal 

Total 


Material 

0. 010 in. 

0. 020 in. 

Components 
Independent 
of Sheet 
Metal 
Thickness 

CR Steel 

i 


1. 13 

CR Steel 

i 


0. 47 

Stainless Steel 

0. 71 

1. 42 


Stainless Steel 

0. 57 

i. 14 


CR Steel 



0. 42 

TZM 

Molybdenum 



0. 71 

TZM 

Molybdenum 



1. 96 

CR Steel 

0. 51 

1. 02 


CR Steel 

0. 08 

0. 15 


Alnico V 



1. 50 




0. 49 


1. 87 

3. 73 

6. 68 


8. 55 

10. 41 









Based on the above criteria, the thruster shown in the upper 
portion of Fig. 59 was designed. The nominal current performance is ■ 
1200 mA (beam) at I = 3500 sec. The total input power to the thruster 
as a function of specific impulse is shown in Fig. 60. 

Early in this program an opportunity arose to test a thruster of 
this size and approximate design (constructed with HAC IR&D funds) 
with an oxide cathode of the type proposed for use on this thruster. 

Stable operation was achieved on the first attempt, with efficiencies of 
approximately 200 eV/ion discharge losses at 85% mass utilization. 

This preliminary test verified that the design of the thruster is basically 
sound and suitable for the proposed usage. 

The 30 cm oxide-cathode thruster design for this program is simi- 
lar to the HAC IR&D thruster which uses a divergent magnetic field, 
similar to the SERT II thruster, but maintains the ionization chamber 
length at 15. 2 cm in accordance with the findings of Richley and _ 

Kerslake and HRL Contract NAS 3-9703. Other modifications are the 
incorporation of a center support between the accel and the screen 
electrodes, and accel electrode insulator supports attached through a 
double leaf spring arrangement. This support was designed to absorb 
the bending torque which would be transmitted to' the screen electrode 
as a result of the differential in thermal expansion. The propellant 
enters the ionization chamber through'two vapor lines, which feed into 
a distribution plenum located within the chamber. The propellant passes 
into the ionization chamber through twelve outlets in the plenum which 
are covered with a fine mesh screen. This configuration reduced the 
thruster over-ail length and the possibility of external leakage of the 
mercury vapor. Leaks which might occur between the plenum and the 
back plate remain within the ionization chamber. 

The thruster structure has been designed around simple sheet 
metal fabricating techniques' irf order to reduce thruster fabrication costs. 
Thus, some of the complex structural configurations which might be 
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formed by spinning were formed from sheet metal and fusion heliarc 
welded to’ form the desired shapes* For example, the propellant mani- 
fold mentioned above could be formed by spinning or drawing a single 
sheet* This would require expensive tooling. Instead, this part, which 
can be loosely described as a double flanged annular ring with a cross 
section shaped like a squared off capital Greek letter omega, was 
formed by rolling two different size hoops from strips and adding three 
rings as shown in Fig. 61. The tooling required for this operation 
were simple lathe- turned rings and were used to hold the shapes during 
the welding. The method is dependent upon the skill of the welder in 
assembling an acceptable part, but allows changes to be made without 
the loss of expensive tooling. Other structural parts of the thruster 
were fabricated by similar procedures. 

C. Accel Electrode Support 

-A novel accel electrode insulator support configuration was 
designed for this thruster. It is shown in Fig. 62. The insulator is 
completely contained within the cylindrical pieces in the center of the 
structure. The cylinders form shields to prevent the coating of the 
insulator with mercury. The two flat parallel vertical members shown 
to the left of the insulator are leaf springs. Another pair, one of which 
is visible in the figure, is located on the right side. They are attached 
to a yoke arrangement which is shown with four holes located on the 
tang which projects away from the viewer. This tang is normally used 
to attach the insulator assembly to the screen electrode. The stud 
projecting from the top of the insulator shields attaches the accel elec- 
trode to the assembly. During engine operation, differences in thermal 
loading between the screen and the accel electrodes cause a variation in 
thermal expansions. The leaf spring members will flex equally and 
maintain the two interconnected surfaces parallel. A small decrease in 
the spacing of the interconnected surfaces will occur, but these are 
insignificant compared with the changes that occur as a result of thermal 
buckling of electrodes. 
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When the thruster is tested on the thrust stand the accel electrode 
will be mounted on the motion generator. The same insulator support 
will be used, except that it will be inverted from the attitude shown in 
Fig. 62. The tang that attaches to the screen electrode will attach to 
the top rings of the motion generator which encircles the thruster. 

D. Thruster Electrodes 

The thruster electrodes were fabricated by match drilling the 
electrode pair on a numerically controlled drilling machine. The 
orientation of the electrodes is indicated by index marks in the form 
of arrows engraved on the outer edge of each electrode. Match drilling 
of the electrode pair was determined to be the only way that the several 
thousand pairs of holes in the two electrodes could be held in absolute 
alignment with any reasonable degree of assurance. Other methods may 
have held the absolute location of each hole within the desired 0. 003 in. 
nonac cumulative tolerance {which was approached by the tape drilling 
machines used) but would -not yield the results that match drilling 
obtained within a reasonable cost. 

The' electrode mounting geometry is that which has been adopted 
for all of the 30 cm thrusters built to date at the Hughes Research 
Laboratories. The screen electrode is 0. 030 in. thick and the openings 
are 0. 188 in. in diameter. The center- to- center distance is 0. 211 in. 
on a hexagonal pattern, resulting in an open area of 72. 2%,. The accel 
electrode openings are 0. 140 in. in diameter with an open area of 40. 5%. 
The accel electrode is 0. 100 in. thick. The openings in the screen 
electrode are not tapered. 

The material specified for the electrodes is arc-cast low-carbon 
molybdenum. The material as received from the producer was annealed 
at 1650°F {900 C) for 1 hour. A 0. 020 in. bow was observed in the 
screen electrode and a 0. 004 in. bow was observed in the accel electrode 
after machining. The electrodes were fixed between two Blanchard ground 
steel plates for stress relief and flattened at 900°C in a vacuum furnace 
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for 1 hour. The fixture consisted of two heavy (125 lb each) Blanchard 
ground steel plates coated with a 5% weight solution of Baymal in water 
as an antiseizing agent. A 0. 005 in. thick molybdenum sheet coated on 
both surfaces with Baymal was- placed between the plates and the elec- 
trodes and between the electrodes. The electrodes were also coated 
on both sides. The antiseizing agent and the thin molybdenum sheets 
were used to prevent "cold" welding in vacuum. The interface between 
molybdenum pieces probably would not be subject to "cold" welding at 
this temperature, but the steel plates require the antiseizing coating. 

After stress relieving and flattening, the electrodes were flat 
within 0. 002 in. over their entire surface. 

E. Center Support 

A drawing of the center support is shown in Fig. 63. It projects 
into the upstream or ionization chamber side of the electrodes. A 
minimum cone-shaped nut is placed in the downstream or beam side to 
reduce the possibility of beam interception. The cup is attached to the 
screen electrode by a 90° included angle cone shaped shoulder on the 
bushing which screws into the cup, through a matching countersunk 
center hole in the screen electrode. The end of the bushing fits flush 
with the surface of the screen electrode. The insulator fits freely in 
the- cup and bushing and is attached to the bottom of the cup. Threaded 
end fittings were swaged over the ends of the solid alumina rod insulator. 
In addition to being a mount for one end of the insulator, the cup also 
acts to shield the insulator from surface contamination. The other 
threaded end fitting of the insulator passes through a clearance hole in 
the accel electrode and is held by the cone-shaped nut. The metal parts 
of the center support are fabricated of tantalum because of its elevated 
temperature properties and room temperature ductility. 
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Fig. 63. 


Center support. 


F. Accel Electrode Motion Generator and Monitor 

A mechanism is required to move the accel electrode with 
respect to the screen electrode in a known matter to verify the 
resultant change in thrust vector. The mechanism must be able to 
move the accel electrode axially, rotationally, and laterally; it also, 
must be able to cause the accel electrode to tilt with respect to the 
screen electrode. The mechanism must be able to repeatedly produce 
these small motions. (Total excursions : 0. 020 in. axially, axial tilt, 
laterally and rotationally at the outermost opening of the electrode). 

For axial motion it was initially proposed to use three solenoid driven 
stepper motors driving a rotary cam. The accel electrode was to be 
supported through these cams. Rotations of the cams would cause the 
accel electrode to tilt. Several problems rapidly become apparent with 
this method'. ' First, with -the o.v.e rerunning clutch cam drive of the stepper 
motor designed at HRL, each clutch performed, slightly differently from 
the other. The result was that for each pulse to the stepper motor sole- 
noid, a uniform axial displacement of the entire accel electrode was 
not possible. Second, it was not possible to maintain an exact location 
of the accel with respect to the screen, without a very complex electrical 
insulation with unidirectional flexibility and rigidity along the other two 
axes. Ideally, the axial motion generator should function through a 
single stepper motor and be capable of motion in the axial direction and 
not along the other two axes. 

A multiple spring column mounted mechanism was designed as 
shown in Fig. 64. Two plates are employed. They are attached to each 
other through eight flat leaf spring column supports. The leaf springs 
are slightly buckled and held in this buckled attitude by a light wire 
cable, which passes over sheaves located at the center of each of the 
leaf springs. Two lengths of cable are used, each threaded through 
four of the leaf spring sheaves. Two ends on one side are attached to a 
cable end which is fastened to the eccentric* on the end of the stepper 
motor shaft. Of the other two ends of the cable located directly across 
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the plates, one is fixed to the plate. The other end is attached to a 
cable end which is again fastened to an eccentric on the end of the shaft 
of a second stepper motor. Actuation of the first stepper motor in 
effect causes the active length of the two cables to be decreased, 
which in turn pulls in the center of the leaf spring columns. _ This 
causes the height of the columns to decrease in a uniform manner. 
Actuation of second stepper motor causes a buckling of four of the 
columns on the one cable. This causes the tilting of the accel electrode. 

The leaf spring columns used are fabricated from 17-7 PH stain- 
less steel in the spring tempered condition. They are. 0.008 in., 
thick and a nominal 0. 500 in. wide. Eight of the columns are located 
radially around the plates. A 0. 020 in. decrease in column height re- 
quires 2. 2 lb tension in the cable. The rigidity of the leaf spring in 
the width di-r-ection is on the order of three o'rders of magnitude greater. 
Two of the leafs will support 7-1/2 lbs and deflect 0. 001 in. along the 
width direction. 

The stepper motor which was finally selected is a commercially 
available, solenoid operated ratchet drive device. Each actuation 
produces 30° of rotation on the shaft. A minimum 6 in. of torque is 
available at the shaft. For a total excursion of 0. 20 in. or 0. 10 in. 
arm length of the cam, 3. 75 lb or more of force is available. Because 
of the long delivery time of this stepper motor, a motor was designed 
and fabricated using a commercially available solenoid and overrunning 
clutches shown in Fig. 65. The test of this stepper motor at elevated 
temperatures resulted in the rapid failure of the overrunning mechanism. 
The feature of this design Was the capability of adjusting the cam incre- 
mental rotation per pulse. 

The ratchet drive stepper motor had to be modified for this appli- 
cation. The initial assembly as received from the supplier used oil 
filled bronze bushings for the ratchet shaft. These were replaced with 
lightly lubricated instrument type ball bearings. The stop mechanism, 
which consisted of a light flat spring, was replaced with a heavy 
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machined arm that was spring loaded. An additional stop was added to 
the pawl mechanism on the solenoid plunger. This was necessary to 
lock the ratchet when the drive pawl, was actuated and retracted from 
the ratchet. This stepper motoris shown in Fig. 66. 

A third plate is included in motion generator assembly. This 
plate holds the accel electrode supports and is used to supply the 
lateral and rotational movement of the electrode. Two stepper motors 
identical to those used for the axial motion are mounted to the middle 
plate (the upper plate of the axial motion generator). ■ The end of the 
shaft of- the stepper motors is a circular eccentric. The end is attached 
to the third plate through a cam follower which consists of two plates 
with an elongated hole. Two plates are used so that the cam follower 
clearance can be adjusted. The follower for the rotational motion of 
the accel has an elongated hoTe that is parallel to the radius. of..the 
accel.' The lateral motion follower has an elongated hole that is per- 
pendicular to the radius of the acceL Rotation of the cam on the step- 
per motor will cause the cam, to bear against the straight sides of the 
elongated holes and impart motion to the electrode. The elongation 
of the cam followers allows the system to slip when the other is 
actuated. The third ring is held. in place by spring loaded guides 
located along the outer edges and on the top of the plate;, it rests on 
eight teflon disk bearings of 3/16 in. diameter. 

The motion monitor-is a commercially available variable induc- 
tion transformer. The test of the first device obtained has shown that 
the best sensitivity of this device occurs with a 400 Hz primary power. 
The sensitivity obtained was 1. 5 mV/V/0. 001 in. The device is not 
sensitive to radial displacement of the movable core nor to elevated 
temperatures. One monitor will be used to indicate axial position of 
the accel electrode and a second will be positioned to sense its lateral 
and radial displacement. A check of the 400 Hz power available at the 
Research Laboratories indicated that the control of the source was not 
as good as one would desire for this application and also the cost of 
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Fig. 66. Modified solenoid stepper motor. 
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operating this unit for such, a small power requirement, would be 
prohibitive; two alternatives were open: (1) build a 400 Hz power supply, 
which is not a trivial task; or (2) use 60 Hz. The 60 Hz was especially 
appealing after consultation with the manufacturer disclosed that a simi- 
lar transformer, with a higher sensitivity at 60 Hz than the unit tested, 
could be obtained. These have been incorporated for use on the test 
stand. 

G. Electrode Thermal Tests 

Thermal bench tests were performed to simulate the thermal 
load on the electrodes. The first and second tests were performed with 
only-the screen electrodes; in the first, the screen electrode was 
mounted with the edges free except at one point. The thermal gradient 
across the face of the electrode was maintained at- approximately 33°C. 
The second test was performed with the screen electrode and the edges • 
clamped’. The third test was performed with the accel electrode in 
place, the center support attached, and the screen electrode edges free 
except at one point. 

The main mounting of the electrode was the screen collar of the 
thruster structure, which was supported by three machinists angle 
plates positioned on a large granite surface plate. A screen pole piece 
was attached to the collar, and the screen electrode was placed on top 
of it exactly as it would be on the thruster. The assembly was leveled 
by placing a ground gauge block on the screen electrode near each angle 
plate support and checking its height on the alignment scope. A circular 
concentric three element electrical heater was located about 4 in. below 
the screen electrode. A cylindrical radiation shield of approximately, 
the diameter of the anode was placed around the heater and extended 

toward but not touching the electrode. Each of the three elements was 

\ 

separately controlled through a variac. 
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The results of these tests are shown in Figs. 67 through 69* 

These figures show the difference between the heights of the ends of 
screws attached to the screen electrode along three diameters for 
tests 1 through 3. Figure 70 shows the height of a reference block of 
tungsten placed adjacent to the screws protruding through the openings 
in the accel electrode in test 3. The figures show a "potato chip" unstable 
effect of thermal loading on"flat" plates. The electrodes were stress 
relieved and flattened before the tests were set up. In setting up. the 
test the electrodes were- heated and handled and may have been subjected 
to high localized stresses which were sufficient to initiate local distor- 
tions of varying magnitude and direction. 

After the tests were completed the electrodes were inspected 
for any deformation that might have occurred. A deformation of 4 mils 
was measured in the accel and 12 mils in the screen. 
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Fig, 68. Screen electrode thermal deformation. 
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Fig. 70. 


Accel electrode thermal deformation. 
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VII. THRUSTER TEST RESULTS 

The thruster was installed in the 9 ft chamber and tested. The 
first test, performed with eight 13/64 x 25/64 in. rectangular -cross 
section permanent magnets, showed a discharge instability which was 
attributed to a high magnetic field. Higher than normal accelerator 
drains were also observed, but they were later attributed to misalign- 
ment of the electrodes. The engine was removed and two of the eight 
magnets were removed. The alignment of the accel electrode to the 
screen electrode was checked and corrected, and the thruster was 
reinstalled in the chamber. 

Because of the increased spacing. between the magnets, the 
azimuthal variations in magnetic field increased to the point that inten- 
sity variations in the discharge plasma could be observed visually. The 
thruster performance was erratic and relatively poor. 

The permanent magnets were replaced with twelve electromag- 
nets, and the thruster was again tested. With this configuration the 
thruster was operated at 1A of beam current at a propellant utilization 
efficiency of 80% and an ion chamber performance of about 200 eV/ion. 
Ten permanent magnets with a cross section of 9/ 64 by 25/ 64 in. were 
installed, and stable operation at beam currents ranging from 900 to 
1400 mA was measured at accel drains of 0. 4 to 0. 5%. The'discharge 
chamber performance during this test was 220 eV/ion at a mass utiliza- 
tion of 85 ±5%. The accelerator and screen temperatures measured 
during this run at the periphery of the electrodes were 225°C and 
228°C, respectively. 

A second thruster test was performed with the thruster mounted 
within the motion generator structure. The axis of the thruster in this 
case was horizontal, which is not the normal attitude of the thruster on 
the thrust stand. The motion generator caused the accel electrode to 
shift ~0. 002 in. downward. This shift was corrected by replacing the 
spring loaded guide screen with a solid guide screen. 
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The thruster was placed in the chamber and allowed to pump 
down overnight* Application of high voltages caused severe arcing to 
occur between the thruster and the motion generator structural rings* 
Arcing was especially noticeable near the anode insulators, which 
project approximately 0* 5 in. from the surface of the outer shell* An 
electrical continuity test between the motion generator structure and 
ground showed that the structure was electrically floating due to the 
anodized surface of the aluminum structure. The anodizing effectively 
insulated the motion generator structure from the ground. Bonding 
straps were employed throughout the structure to ground and the thruster 
retested. Arcing again occured when high voltage was applied. The 
thruster was left overnight in the chamber to allow it to outgas. Although 
.the frequency of the- arcs appeared to dimmish and operating voltages of 
+3 kV and -2 kV could be applied, stable operation, could not be attained. 
The test was terminated, with the conclusion that the insulating effect of 
the anodizing was causing surface charging and that the arcing was the 
result of the charge buildup. The anodizing was stripped from the sur- 
face of the motion generator structure. Further tests were deferred, 
until the thruster and the thrust stand were* completed in order to 
conserve funds and time. 
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VIII* THRUST STAND DESIGN 


The thrust stand design selected to perform the task of evaluating 
the extensive analysis performed was the floating suspension system 
analyzed earlier in the report (see Fig. 59). Separated into three major 
functional categories, the major component groups are the float, the 
sens or -force motor and the thruster-motion generator. Of the three 
categories, the thruster-motion generator group is the only one that is 
totally floating. The other two groups must, by necessity, be referenced 
to a ground plane. The float is referenced through the pool of mercury in 
which it lies, and the sensor-force motor is referenced through a fixed 
platform. This platform performs the double function of supporting one 
half of the sensor and forcemotor assemblies and provides a horizontal 
reference plane. The design is such that the only physical interconnec- 
tion between the floating and fixed parts of the thrust stand assembly 
could be through the mercury pool and the electrical lines which would 
supply the power to the various floating components of the thrust stand, 
thruster and the motion generator. The mercury pool formed a very 
compliant interconnection to the ground plane; this made it a major task 
in the fabrication of the thrust stand to find an equally compliant and 
flexible electrical lead to go between the ground plane and the floating 
mass. Considerable effort was expended on this task to assure that the 
spring rate for the electrical leads would not be detrimental to the 
performance of the thrust stand. The various components were designed 
and their performance checked to assure that they would be adequate to 
perform the task for which they were intended. The areas of critical 
assembly are defined. 

A. Float 

A flat-bottomed, straight- sided annular float is used to support 
the floating mass of the thrust stand on a mercury pool. The configura- 
tions of the float and of the pan were chosen to allow a central passage 
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for the electrical leads and to provide a well for the thruster feed system. 
The volume of the float is such that it will support 28. 8 lb/in. of immer- 
sion in the mercury. The clearance between the float and the pan walls 
is 1-3/8 in. at the surface of the mercury. This was necessary to 
minimize surface tension effects by providing separation of the mercury- 
metal interfaces should mercury adhere to the walls of the float or pan. 
Approximately 5/16 in. below the surface of the mercury, the clearance 
between the pan and the float is reduced to 1/2 in. This reduces the 
volume of mercury required without affecting the sensitivity of the 
thrust stand. The sides of the float have been lined with teflon to in- 
hibit the adherence of the mercury to the walls of the stainless steel 
float. < 

B. Electrical Conductor Cables 

The only mechanical link between the floating and the fixed mass 
of the thrust stand are the electrical conductor cables required to trans- 
mit electrical power to the various parts of the thrust stand and the 
thruster. A minimum of forty separate conductors, of varying current 
capacity, will be required. The electrical cable with the most flexibility 
and best current carrying capacity was made by Cicoil Corporation. Of 
the various samples supplied, the most promising appeared to be an 
eight #32 AWG conductor flat cable with a silicon rubber insulation. The 
conductors were specially stranded with 64 strands of #50 AWG 0. 001 in. 
wire. Eight 8-conductor cables are used with a total of 64 leads #32 AWG 
size. Where the current exceeds the current carrying capacity of the 
small wire size, several conductors are used in parallel to increase the 
current carrying capacity. 

A test fixture was fabricated for measuring the low forces, such 
as the resistance to flexing of the electrical cables between the fixed and 
the floating portion of the thrust stand. The test fixture shown schemati- 
cally in Fig. 71 is a one degree of freedom air bearing consisting of a 
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Fig. 71. Air bearing, cable spring rate test fixture. 


triangular cross section plenum perforated on two surfaces with, many 
small holes and a slider which fits closely over the perforated surfaces. 
Pressurizing the plenum with air or nitrogen causes the slide to lift, 
and a virtually friction free system results. Tilting the plenum will 
cause the slide to move to the lower end of the plenum. With no 
restraints the amount of tilt required to start the slide moving is less 
than 0. 000 5 in. on the 12 in. long plenum. The motion of the slide is 
detected by the optical sensor designed for the thrust stand. With the 
sensor electronic circuitry employed during this test, a sensitivity of 
0. 0018 in. per 100 mA is obtained. With a slide weight of 250 g, 
including the sensor aperture, counterweights, and half of the cable 
weight, it is possible to determine the spring rate of the cable by 
measuring the displacement of the slide as the air bearing plenum is 
-tilted .through a few degrees of arc from the following: 

, W • Ah 
k = L * i ‘ 454 • 

where 

k = spring rate of the cable, in pounds per inch 

W = weight of the slide, -in grams 

Ah' = displacement of- the airbearing plenum on one 
end, in inches 

15 = length of the plenum between pivot and microm- 
eter tilt mechanism, in inches 

t, s displacement of the slide from null to equilibrium 
with W ♦. Ah/Li force. 

A test of two cables with opposing 90° twists indicates that this 
cable configuration results in a spring rate of 0. 0 5 lb/in. Because the 
thrust stand at null will, remain well within 0. 0005 in. of its initial 
position, the forces from the leads will be 

F < (5 x 10“ 2 lb/in. ) (5 x 10" 4 in. ) (8 leads)' = 2 x-10' 4 lb, 
which is comparable to the force to be measured. 
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Several conductor cable configurations ihave been tested; with 
the results shown in Table XIV. Because the cable 'must be able to move 
freely in. all three axes, the looped cable is the most desirable configu- 
ration. The direction sensitivity shown between tests 8c and 9 is unde- 
sirable,. It was thought that- splitting the insulation between conductors, 
as in test 4, would reduce the effective "web" thickness and result in a 
direction insensitive flexible cable; however, the test results do not 
support this assumption. These show instead an inconsistent spring 
rate which is attributed to the random rubbing of the strands against 
each other. 

A helically wound cable configuration would result in a low spring 

rate and be insensitive to direction. The supplier of the cable was not 

able to provide specific instructions for coiling or forming the cables. 

Tong duration immersion of coiled cable in solvents (such as acetone, 

alcohol, or trichlorethylene) or heating at 260°C were suggested as 

possible means of setting the cable. Further -investigation suggested 

the possibility of chemical catalyst and sulfur vulcanization at the 

7 8 

temperatures of the previous air firing. ’ The first method tried was 
to immerse the cable in methyl- ethyl- ketone (MEK) after heating at 
280°C for 30 min and then reheat to 280°C to drive out the MEK. In 
one case -the coiled cable had a small spring-back when released from 
the mandrel. A second trial resulted in blistering of the insulation. 

Sulfur vulcanization by sprinkling flowers of sulfur directly on the 
coiled cable caused the cable to adhere to the mandrel. 

Heat aging of the cable at 260°C in air yielded the most consistent 
results. Several samples were prepared and aged for varying lengths of 
time. The best results were obtained after a 4 hour aging period; shorter 
aging resulted in partially set coils, and longer aging caused the insulation 
to adhere to the mandrel and to soften and tear readily. 

The 4 hour sample retained its resilience while holding. the coil 
"set. " It could be stretched -50 to 60% and bent 45 to 60° without losing 
its set; stretching or bending beyond these limits causes deformation. 
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Tests of the helically coiled cables with the coil ends clamped 
perpendicular to the wires show that the spring rate of this configuration 
-is not direction sensitive and is 0. 051 lb/in. , which is slightly greater 
than the spring rate of the cable in the direction perpendicular to its 
width. 

Close inspection of the coiled cable movement while it was 
clamped as described above revealed that although five turns of cable 
were held in the space between the two clamps, most (70 to 80%) of 
the flexing occurred within the turns adjacent to the clamps. A clamp- 
ing fixture was prepared for holding at least one turn of the cable. In 
this fixture 2-1/2 turns of the helically formed cable are held against 
the walls of a cylindrical cavity by means of an expanding cylindrical 
spring, inside the helix. This device imparts a stiffness to the-fi-r-st 
free turn, which would cause all of the turns in the cable to be flexed 
equally and result in a lower spring rate. It was found that this method 
of holding. the ends of the cable reduced the spring rate of the cable in 
the 0°, 45°, and the 90° directions to 0. 003 to 0. 002 Ib/in. This is a 
reduction of one order of magnitude over the previously reported value 
•for the helically wound cable with the clamps holding the cable ends 
perpendicular to the direction of the wires. 

For eight cables and a sensor sensitivity of 0. 0005 in. , the 

0. 003 Ib/in. spring rate would result in a resisting force of 1, 2 x 10 ^ lb, 

which is one order of magnitude less than the measured minimum force 
-4 

component (10 lb). 

Dielectric 'tests of the formed conductor cables were made to 
determine their capability at high voltage. High voltage applied on the 
outermost lead resulted in an occasional breakdown to the surrounding 
cylindrical holder. This would occur at a voltage of around 1 kV. Insu- 
lation breakdown between leads would occur at a voltage greater than 
4 kV. The cable clamping holes in the fixed and floating platforms were 
sleeved with ceramic tubes which projected 1/8 in. above the surface of 
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the platforms. The sleeves were epoxied into place_ and provided the 
additional insulation required to insure that arcing would not take place 
at these points. The spring clips within the cables were wired to the 
terminals to avoid an electrically floating piece at that point. 

Four of eight cables were formed with clockwise helices; the 
remaining four were formed with counter-clockwise helices. With one 
or two exceptions the circuit was completed through the same cable to 
reduce the magnetic interactions which might result. 

A. wiring schematic showing wiring arrangement at the floating 
base is shown in Fig. 72. 

C. Force Motors 

The force motor dfesign goal was 1/2 lb of force for an input of 
1/4 A at 15 V. The thrust stand analysis indicated that a force of 0. 11b 
would be adequate, and that the higher force capability of the force 
motor would enhance the response of the system. The force motor 
design is, similar to that of a voice coil of a loudspeaker, except for 
the air gap. In a normal loudspeaker the air gap is very small and 
the wire wound coil is usually mounted on a thin backing of paper. The 
coil is then centered on the airgap by a diaphragm or a flat spring. 

The force motor must be capable of accommodating both axial and 
lateral displacement, and thus by necessity has a larger air gap than 
the loudspeaker voice coil. The effect of this large air gap is then a 
significant increase in the flux leakage at the air gap. 

A design flux density of 2000 G was assumed for a 0. 25 in. air 
gap of the magnetic circuit. A coil diameter of 1. 50 in. was selected 
based on the flux of Alnico V and the size of readily available plug mag- 
nets. The coil width of 0. 375 in. was selected to overlap the pole faces 
of 0. 313 in. by 10% on each side. The 0. 313 in. pole face width at 
1. 438 in. diameter matches the cross sectional area of the plug magnet. 
Similar area matches were made throughout the iron. The magnet 
assembly is 2-1/2 in. in diameter and 1-1/2 in. high. 
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Fig. 72. Thrust stand wiring. 
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Initial tests of the force motor showed a force of 0. 13 lb at 
300 mA, By magnetizing the magnets in situ, a nearly four-fold increase 
in the force, to 0, 5 -lb at 300 mA, was realized. 

r r 

A z axis force motor assembly is shown in Fig. 73. The x-y 
force motors are identical to. that shown, except for the mounting. The 
z force motors are mounted with their axes vertical and the line of force 
is in the vertical direction. The x-y foi;ce motors are mounted with 
their axes horizontal; the axes of the two x force motors are perpendicular 
to those of the two y force motors. 

A guide ring is pressed on the body of the magnet assembly. The 
mount of the coil assembly has a machined shoulder which overhangs the 
coil and fits within the guide ring of the magnet. The clearance between 
the diameter of the shoulder and the guide ring is 0. 031 in. and that 
between the coil and the magnet pole faces is 0..040 in. Thus the guide- 
shoulder acts to protect the coil by acting as a stop. It is also intended 
as an aid in the installation of the force motor during thrust stand 
assembly. The force motors must be positioned on the tables as an 
assembly. If a 0. 03 in. thick split ring is placed in the clearance 
between the shoulder and guide ring and 0. 031 in. pins are placed through 
four holes in the guide rings to maintain the spacing between the bottom 
of the shoulder and the top of the magnet subassembly, the force motor 
assembly is in the null position. The split ring and spacer pins should 
be left in place until all of the force motors are in position. 

The curve in Fig. 74 is representative of force motor performance. 
These data were obtained by placing. the sensor, force motor, and elec- 
tronics from loop #6, one of the z-axis loops, on the test fixture shown 
on Fig. 73. 

D. Displacement Sensors 

A displacement sensor which utilizes a light beam and a pair of 
matched photocells has been developed (see Fig. 75); the light falling 


185 



EI080-I3 



Fig. 73. Axial force motor assembly. 
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Fig. 75. Optical linear displacement model. 
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onto the cell is apertured by a slit mounted on the floating table. The 
relative motion between the slit and the light bulb and sensors which 
are mounted on the fixed table causes a variation in the amount of 
light reaching each half of the split photocell. This produces a difference 
in voltage across the matched load resistors which may be interpreted 
directly in terms of the relative motion. 

A system of this type using a United Detector Technology photo - 
detector and gallium arsenide light source has been designed, fabricated, 
and tested. As shown in Fig. 76, this model, actuated by a micrometer 
drive, has demonstrated resolution of linear motion of less than 0. 0001 in. 
and good sensitivity. This high sensitivity and the absence of any me- 
chanical coupling between the floating and fixed suspension make this a 
very attractive system for this application. 

E. Count e r weight 

In the translational mode of accel electrode-misalignment the 
accel electrode will be shifted off center in incremental amounts, caus- 
ing a corresponding small shift in the location of the center of gravity 
in the horizontal plane. This shift in the horizontal location of the c. g. 
will be compensated by the counterbalancing mechanism. The mechanism 
consists of two motor- driven weights which move along the radials of two 
of the four z force motors. These radials are at right angles to each 
other. The motors are reversible ac. instrument type motors. The 
output shaft rotates at ’6 0 rpm. The ten tooth gear drives a 30 tooth 
gear on the end of a screw shaft. The shaft then rotates at 20 rpm and 
the screw threads advance or retract a 3/8 lb weight, dependent upon 
the direction of rotation of the screw shaft. The. 1/32 in. pitch thread 
on the screw shaft causes the weight to move 0; 031 in. per revolution. 
Thus when activated for 2 min. , the weight will move 0. 625 in. and a 
change of 0. 234 in. -lb can be made. The counterweight mechanism is 
shown in Fig. 77. To prevent damage of the gears or of the motor, limit 
switches are provided at both ends of the weight travel. The total travel 
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Fig. 76 


Sensitivity of linear displacement sensor. 








of the weight is 1 in. The response of the motor to start and stop is 
within a few milliseconds and thus with the low speed of the weight 
relatively small adjustments can be made. 

F. Sensor-ForceMotor Electronics 

The proposed electronic control loop configurations were dis- 
cussed previously, and a detailed analysis (paper study and analog 
simulation) of the expected system performance was presented. 

The electronics for the single control loop, which was con- 
structed and tested, is shown in Fig. 78. The component values and 
the amplifiers indicated are those suggested (as nominal) for each of 
the eight control loops to be used in the final system. 

Although no completely realistic model was available of the 
actual dynamic load which the- force motor will drive in the final sys- 
tem, various closed loop tests were performed using the electronics 
of Fig. 78 and the mechanical setup shown in Fig. 79. The tests were 
intended to determine qualitatively the compatibility of components, the 
ease of compensation adjustments and their relative effects, the pro- 
cedure for simultaneous mechanical and electrical adjustments, etc. 

The adjustable compensation gain constants functioned as 

expected. The loop was unstable with only position information and 

became stable as the derivative (rate) compensation was increased. 

Error displacement could be driven to essentially zero (zero average 

value) by providing a signal which is proportional to the integral of the 

error. With all gain constants (K , K,, KJ adjusted to apparently 
* j P 

satisfactory values for reasonably well damped transient response, no 
electrical z'ero drift was detected over a period of appro xim ately 6 
hours at ambient room temperature. 

The indicated displacement noise (due to vibration of ’the build- 
ing, etc. ) was approximately 30 MV peak-to-peak (at the output of the 
sensor electronics) and corresponds to a displacement of ±25 pin. It 
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is of more importance that in this particular test setup the background 
vibration produced a force motor current variation equivalent to ap- 
proximately ±5 mlb. The wave shape at the sensor tended to be sinus- 
oidal (in the frequency range of 30 to 60 Hz) and the amplitude was 
dependent on the total weight placed on the balance, decreasing with an 
increase in weight. It should be emphasized that this sinusoidal dis- 
placement had very little effect on the indicated differential force when 

read on a current meter with slow response. A static differential force 
-4 

in the range of 10 lb could be accurately determined. Based on the 
above and other pertinent information on the expected vibrational ef- 
fects, a low pass filter will be incorporated in the current/force 
ins tr umentation. 

G. "Vibration 

With the aid of a very simplified model, the relationship between 
the vibrational displacement of the base of the thrust stand and the 
restoring force motor drive current will be discussed. 

The simplified model to be used is shown in Fig. 80. In this 
single degree of translation model, Xj represents the position of the 
outer housing of the thrust stand (termed the pan) located on a vibrat- 
ing structure and yL^ is the relative position (with respect to a fixed 
reference frame) of the floating portion of the thrust stand (called the 
table). The damping and spring coefficients are denoted by B and K, 
respectively, and F.^ is the restoring control force which is of course 
proportional to the difference between Xj and X^. F^, is the thruster 

force component associated with the model and is to be determined 
for a knowledge (measurement) of F^. M denotes the total mass of 
the table and thruster. Although the model is highly idealized, it is 
sufficiently realistic to give considerable physical insight into the 
vibration problem and indicates possible solutions. 
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It can be assumed that Xj is a known quantity since it can be 
derived from experimental measurements of the acceleration spectrum. 
In an attempt to follow (with some simplification) the notation used 
previously, let 


KM ? 

F e (S) = (S* + aj S + a 2 ) 


G U (S) 


1/M 


B_ q . K 
M S + M 


where 


a l 

a 2 = Kj/Kjj, rad/sec^ 

P = force motor gain, lb/V 

6E x 

— — = displacement transducer gain, V/in. 
°x 

K = proportional gain, V/V 
P 

= derivative gain, V-sec/V 
Kj = integral gain, V/V-sec. 


U rrf ) k d« rad/sec 


K /K~, rad/sec 
p D 


From Fig. 81 and the above notation, the block diagram shown 
in Fig. 82 can be developed. For = 0 the configuration reduces to 
one dis.cussed previously which assumes no external disturbance force. 
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The effects of can be viewed equivalently as a disturbance force 
term which is superimposed on the thruster component F^. for a thrust 
stand in a vibration free environment. Therefore, Fig. 82 is an 
equivalent representation of the block diagram of Fig. 81 where the 
disturbance force Fy is given by 

F y (S) = [SB + K + G e (S)] Xj(S )‘ 

= Gy(S) X^S) 

Now assume that the compensation'gains are adjusted such that 

a, = K /K~ = B/M 

1 P JJ 

a 2 = Kj/K d = K/M 

Then G y (S) = (K p + K^S + K^./S) where 

K = K -B + K 
p o 

1L = K M + B 
D o 

K T = K K ; 

I o 

hence Fy is the weighted sum of the integral, the derivative, and the 
scaled value of Xj. 

The closed loop transfer function becomes 

G ( S ) = S + K = St + 1 ’ T = KT 
o o 
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Hence 


F (S) F (S) 

F R< S > = Sr + 1 + St + .1 - 

= f re^ + f rv^ s * 

If Z denotes the- actual vibrational displacement of the support structure, 
then Xj can be given by Z(S) = A(S) X^(S) where A(S) represents the 
mechanical transfer function of the coupling to the thrust stand pan. 
is therefore 


F RV^ 


(K p + K d S + Kj/S) A(S) 

(St + 1) 


Z(S) 


For the chosen model and values of gain coefficients the prob- 
lem becomes one of determining F^^ from a measurement of F^. The 
most attractive methods for obtaining accurate measurements of 
are to reduce by isolating the thrust stand from the vibrating sup- 
port structure (make |A(S)| « 1 for all frequencies of interest) or 

to separate F-^g from F^ by sufficient low pass filtering, or a com- 
bination of both of the above methods, 

A vibration survey of the 9 ft vacuum chamber was performed. 
Various modes of chamber operation were investigated. Significant 
increases in vibration levels were observed during the operation of the 
holding pump for the chamber diffusion pump, at the initiation of the 
filling of the L* cryowall, and the filling of the diffusion pump*LN 2 
baffle. Fortunately all three of these functions can be turned off manually 
during the thrust stand test and thus do not affect chamber vacuum per- 
formance, The significant result of the vibration survey was the 

_ 3 

observation of an acceleration peak of 3 x 10 g occurring at 25 Hz. 
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The frequency range is such that, it may affect the thrust stand operation, 
imposing a background noise to the sensor which cannot be filtered out 
without affecting system response and sensitivity. 

Vibration isolators in the form of air inflated cushions have 
resonant frequencies on the order of 1. 5 Hz which could be employed to 
isolate the thrust stand platform from the chamber. However, because 
there is some concern about employing such a system within the vacuum 
chamber, other alternatives were considered. The best plan, which was 
adopted and accomplished with company funds, was the upending of a 
normally horizontal 4 ft diameter chamber and mounting the entire 
chamber on the air cushion vibration isolator as shown in Fig. 82. The 
chamber was relocated next to the 9 ft chamber because of the higher 
-ceiling clearances in that area. All of the vacuum pump lines, the liquid 
nitrogen lines, and cooling water lines were .equipped with vibration 
isolation to avoid transmission of pump and external exciting forces. 

The general configuration of the chamber is a cylindrical chamber with 
its axis vertically supported 4 ft off the floor and on four columns with 
three vibration isolators at the interface between the chamber and the 
columns. Varying the pressures in the isolators permits small adjust- 
ments in the attitude of the chamber. The access to the chamber is 
through the two end bells. A cable device will lower or raise the end 
bells and will also be useful in raising test items to the upper end of 
the chamber whenever the chamber is used in that mode. 

The thrust stand is mounted on the lower bell on a platform which 
raises it approximately 1 ft above the inside surface of the end bell. 

This allows the electrical connections to be located beneath the thrust 
stand. The platform on which the thrust stand rests can be levelled. 

It will also contain cooling coils which will allow a control of the 
temperature of the mercury pool. The installation is shown in Fig. 84. 
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Fig. 83. Preliminary sketch of the 4 ft diameter chamber in the 
upright position. 
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(a) 


(b) 


Fig. 84. (a) Console with thrust stand in the vertical vacuum chamber, 

(b) Close-up of thrust stand. 
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H. Mercury Level Control 

A direct finding of the integration test of the thrust stand system 
(described later) was the need for some method of controlling the 
mercury level in the pan of the thrust stand or of the location of the 
axial c. g. , or both. The -level of the mercury within the pan con- 
trolled the alignment of the horizontal slot in the sensor for the z axis 
force motors. The location of the axial c. g. should coincide with the 
metacentric height of the -floating mass. For the short term tests the 
amount of mercury used from the feed system will be quite small and 
the change -in the axial center of gravity will be quite small also. The 
change in the metacentric height will be even less; therefore, to com- 
plete this phase of this program it will be necessary only to be able to 
adjust the height of the mercury level to compensate for the change in 
the weight of the system. A bellows mechanism has been designed for 
this purpose. In application, the bellows and the tubing, leading up to 
the mercury pool will be completely filled; as the need for a change in 
the mercury level is dictated by the shift in the null position of the z 
sensors, mercury can be added or removed by changing the length of 
the bellows. This is a manual operation and the bellows will be 
located outside the chamber. 

I. Mechanical Assembly of the Thrust Stand 

The mechanical assembly of the thrust stand was one of the 
major efforts expended in this program. It was not possible nor 
desired to build each component with close tolerances such that the 
components would drop into place. The -parts, in general, were built 
to normal tolerances and suitable adjustments and clearances were 
provided so that critical locations, clearances, alignments, and adjust- 
ments were made during the assembly process with procedures using 
simple assembly aids, shims, fixtures, and precision inspection tools. 
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The force motors are a prime example of the techniques 
employed. The location of the individual force motors with respect to 
each other was critical, but close tolerances were not machined into 
the mounting bracketry or into the large plates to which they were 
attached. Instead, the magnet assembly was carefully/ located on one 
of the' large plates, the fixed plate, of the thrust stand. The angular 
location of the magnet assembly was determined by using a^precision 
rotary table and an alignment telescope to within 5 min. The radial 
location of the magnets was set to. within 0. 002 in. with a micrometer. 
The magnet assemblies were then do welled to the fixed plate. The coil 
assembly was then positioned on the magnet assembly with shims and 
pins as described previously. The coil assemblies were then attached 
to the floating plate of the thrust stand. This in itself was not a trivial 
process. The distance between the floating plate and the fixed plate was 
set and held at 3 in. by ground and hardened spacer blocks. Each of the 
eight force -motor assemblies was held in respective "null" position sub- 
ject to machining tolerances, which must be either lapped or shimmed. 
After all eight of the force motor heights were fixed, the coil assem- 
blies were do welled to the floating plate. The previously mentioned 
gage blocks were left in place during the rest of the assembly 
procedure. 

A similar procedure was employed in positioning the x-y sensor 
slits along the edge of the floating plate. The slits of the x-y sensors 

* i 

are- vertical; thus the mating light source and photocell mounted on the 
fixed plate were easily adjusted. For the x or horizontal slit such an 
adjustment provision could not easily be provided; its adjustment and. 
location had to be accomplished with the light source and photocell. 

Adjusting the axial location of the center of gravity was accom- 
plished by placing the entire floating mass in a horizontal position over 
knife edge pivots and moving the thruster and motion generator- assembly 
or the float assembly. The center of gravity was adjusted to coincide 
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with the. calculated metacenter of the floating mass. As reported earlier, 
such a condition will result in an entirely neutral system with no forces 
that will tend to right it or upset it. This technique- was used instead of 
the calculated center of gravity method because of its greater simplicity 
and accuracy. 

J. Integration Test 

Test results, indicate that the concepts on which the thrust stand 
is based are valid and that the mechanical and electromechanical compo- 
nents and the electronic control circuitry chosen for the multiloop con- 
figuration are -generally compatible. 

Minor electronic modifications were required to achieve satis- 
factory dynamic stability. This was due partially to noisy signal lines, 
variations in the light sources, insufficient power supply and signal 
ground decoupling. In addition, a two-fold increase in the floating mass 
combined with an apparent low sensor gain {compared with early tests) 
led to a loop gain approximately one order of magnitude lower (for most 
loops) than that considered in the initial analysis. The gain loss has 
been partially offset by increasing the sensor light source bias current 
for greater output. 

The electronic removal of an environmental noise (electronic 
and vibrational) was relatively straightforward. First-order instru- 
mentation filters with time constants, of 1 to 2 sec have provided 10:1 

-4 

signal-to-noise ratios when recording static forces in the 10 lb range, 
and can, of course, be improved upon with an increase in instrumentation 
response time. At present it appears that the limiting factor in the sys- 
tems sensitivity may be the drift in electronic null settings. 

Tests of the thruster operating at 0.75 A beam current with 
over-all beam voltage of 4 kV (2.4 kV screen and neg 1. 6 kV accel-) 
have been conducted with the thruster beam directed upward within the 
vertical 4 ft diameter chamber. These tests demonstrated that the thruster, 
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power conditioning, thrust stand and its associated data recording ap-' 
paratus were all functional. The system was operated for several hours 
on all sensitivity range's. The effect of small arcs in the accelerator 
region were observed as abrupt transients followed by a damped oscil- • 
lation of the freely suspended thrust stand. The initial offset of the 
thrust due to manufacturing misalignment was observed. Unfortunately, 
a failure in the feed system made it impossible to gather additional and 
more quantitative data. Because of the delicate nature of the instrumen- 
tation, time was not available to make the necessary repairs and continue 
system testing during this phase of the program. 
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APPENDIX I - DETERMINATION OF GEOMETRIC (VIEW) FACTORS 


A. Definition of Geometric (View) Factors for Radiative Heat 

T ransfer 


£ 992-1 



The definition of the geometric 
(view) factor F^-A^ which 
enters into the calculation of 
the radiant heat transfer between 
two finite surfaces A^, A^ 
follows from the assumptions 
that: 


1. 


2 . 


i 


The intensity of heat radiated from an infinitesimal area 
element dA^ is proportional to dA^ and varies with 
direction as the cosine of the angle from the normal n^ 
to dA i . 

The quantity of heat incident on an element dA^ from 
dAj is proportional to the solid angle dco = dA 2 cos ^2^ T 
subtended by dA 2 as seen from dAj. 


The ratio 


/ / cos Pj cos p z 


A 1 A 2 


dA 1 dA 2 


/ / cos p^ dAj dto 

A , 2ir 


’ A, / 1 


COS COS P 2 


1 A 1 A z 


irr 


dA 1 dA 2 


a F 


A 1 - A Z ’ 
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of the heat radiated by area 
radiated by -A-^ in all directions 
factor from area A^ 

From this definition of F 
relation 


which is incident on area A^ to the heat 
( 2ir solid angle) is defined as the view 
to A^. 

A . , the very useful reciprocity 

A r A 2 • / * 


A 1 Fa ,- a 2 = A 2 F A 2 - A ! 


follows immediately Another useful relation for deriving one view factor 
from others is based on conservation of energy. This states that if a sur- 
face radiates to surface A^ and A^ = A^ + A^, then 


r* — tt i + F n 

A 1~ A 2 A l" A 2 A l' A 2 


B. Four Elementary View Factors 

Formulas for the four basic view factors which are used in the 
determination of the over-all grid system view factors are given below. 


1 . Cylinder Walls to Disk 


In the previous quarterly report the view factor from the 
walls of a cylinder to a disk was derived. With new notation correspond- 
ing to the figure below, the formula is 



2. Section of the Cylinder Wall to Disk 


Using the reciprocity law in the relation. 
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-1 


gives 


F D-2 ” F D-A + F D-1 


A 2 F 2-D “ A A F A-D + A 1 F 1-D* 


where 
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A 1 

= 2iri(H-b) 

A 2 

= 2irl H 
2 


a d 

= irr 

/H-1 

f i-d 

= f wd 

(-T 

/ H 

F 2-D 

— F 

WD 

r- 


The view factor from the section of a cylinder A’ to the disk D is then 
found to be 


F AD (H,b,J,r) = 


H F - - 

b ID 


(H r \ H-b ( H-b r\ 

5 \i » i/ " b WD Vi ’ i/ ' 


The next two view factors are found, in many texts. See, for 

't' 

example, Eckert. 



‘E. Eckert, Heat and Mass Transfer (McGraw-Hill, New York, 1959). 
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Area Element to Disk 



5 . 


Disk to Off-Axis Disk 

E992-6 




The view factor from a disk A. 


A " is found from the integral 


to a parallel but off-axis disk 

A 2 1 

of F dA r A 2 : 

/ 


A r A z 


j. 

A 


_A *^1 

1 k dA l 2 


This integral must be calculated 
numerically; the most convenient 
form is 


a +r, 
o 1 


1 


r 


<1 A A 





S(a) = 2a- cos 


D(a) 


D(a) = 


2 , 2 2 
a + a - r. 
o 1 

2a 


In the applications below, the view factor is referred to as 

F SD (h ’ r l’ r 2’ a o ) ’ i<e " 


a +r, 
o 1 


F SD (h, Tj, r 2> a Q ) a 


i- f 

L 1 J 


a -r 1 
o 1 


F dA r A 2 S ' a > da 


c. Calculation of View Factors for the Grid System 


1. View Factor Fj_^ 

The figure at right 
shows a pair of aligned apertures 
in the perforated plates S and A, 
Node 1 corresponds to the disk 
which is at the back of the plate S 
aperture and node 2 corresponds 
to the cylindrical aperture wall. 
(For simplicity, the screen aper- 
tures are not countersunk. ) The 
view factor Fj_ 2 is given by 


E992 

PH 



1-2 


= 1-F 


1-9 


>- F DD |r ,'V l I l 


FRONT 
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2 . 


View Factor F 


1-7 

Node 7 is a disk which is the front of the plate A aper- 
ture. Auxiliary nodes 8, 10, 11, and 12 are also used in the calcula- 
tion of F^ Nodes 8 and 10 are disks concentric to 11 and 7, and 
nodes 11 and 12 are the remaining annuli. An "aperture" effect is 
present in the calculation of F^ ^ because of the nonzero thickness d^ 
of plate A. Thus F^_^ * s a ^ unc tio n of d^. It would be difficult to 
account for this partial shielding between nodes 1 and 7 in the view 
factor calculation. Instead, upper and lower bounds are found for 
Fj y(d^). From the geometry it is clear that for d^ > 0, 

I 

F l-7< d A> < • 

A good lower bound for F^_^(d^) is found as follows.- Using the equations 

F l-7 (d A ) = F 1 - 1 2 + F l-10 


1-12 


A. . A_ 

— — F + — — F 

A, 11-12 + A * 


1 


1 


8-12 ’ 


we obtain 


F l-10 + Aj F 8-12 + Aj F 1 l-12^ d A^ = F l-7^ d A^ < F l-7^ 


= F 


■^1 1 

+ ^ F, 




1-10 ‘ Aj *11-12^ + A F 8-12 * 


Thus the inequalities 


F l-10 + A 1 F 8-12 < F l-7^ d A^ < F l-7^°^ 


( 1 - 1 ) 
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bound F^_^(d^) below by taking p_ j 2^ d A^ = ® and above by having 

F ll-12^ d A^ = F ll-12^' Since F ll-12^ d A^ < F ll-12^ <<: ** the 
complicated view factor F^ ^(d^) * s ac t ua -Hy negligible. In terms of 

the elementary view factors given, eq. (1-1) becomes 


/r \ 2 

F DD^ r S’ r 10’ + (r^j { F DD^ r A J r A* " F DD* r A’ r l 0’ W * } 


F 1 - 7 ^ d A ^ > F DD^ r S’ r A’ w ^ 


<F l-7< d A> 

( 1 - 2 ) 

(1-3) 


where 

i 


w = dg + d + d^ 


10 


" r A " d+d. 




3. View Factor Fj_^ 

Node 4 corresponds to all of the rraterial on the back 
side of plate A. The auxiliary nodes 9 and 13 are disks corresponding 
to the front of the plate S aperture and the back of the* plate A aperture. 
The view factor F^ ^ will be found from the relation 


F l-4 = F 1 -9 " F 1 — 1 3 1 * 

where node 13’ corresponds to the disks at the back of all of the plat^ 
A apertures. The view factors F^_^ and ^ , are given by 

F l-9 = F DD^ d S’ r S' r S^ 

F 1 - 1 3 1 = F DD^ d S + d ’ r S’ r A^ + 6F SD^ d S + d> r S’ r A’ S o^ * 

In the expression for F^ ^ i only the six apertures surrounding node 
13 were considered, and the shielding due to the thickness of plate S 
was neglected. The approximate expression for F^_^ is then 
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1-4 


= F 


DD (d S’ r S* r S } " ^ F DD {d S + d ’ r S’ r A> + 6F SD (d S + d> r S’ r A' S o ) > * 


4. 

in plate A. 


, View Factor F ^_5 

Node 5 represents the cylindrical wall of the aperture 
It is easily seen that 


F 


1-5 


F 1 — 1 3 " F l-7 


F 


1-5 


I 

where the value 


of 


F 


= F DD* d S + d> r S’ r A^ " F l-7^ d A^ 

(d. ) is chosen from (1-2) and (1-3). 
1-7 A 


5. View Factor F ^_7 
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As in the case of F^_ 7> there 
is partial shielding in the path 
from node 2 to node 7, because 
of the nonzero thickness d^ of 
plate A. Instead of calculating 
F_ (d A ) precisely, an estimate 
is made based on upper and lower 
bounds. An upper bound which 
neglects the shielding completely 
is F 2 _ ? (d A ) < F 2 _ ? (0). A lower 
bound which omits the difficult 

view factor -^ 2-15 ^ 2 _ 7 (d A ) > 

F, ... In terms of the elementary 
2-14 

view factor formulas these bounds 
become 


F A n (w ’ d S’ r 14’ r S ) < F Z-7 (d A > < F AD (w> d S’ r A’ r S } 
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where 


w 


d s + d + d A , 


6 . 


r 14 “ r A d * r S ' r A } * 

View Factor F 2 _ 5 

The view factor F-, c is found from the relation 

L>~ D 


F 2-5 “ F 2-13 F 2-7 ; 

node 13 is shown in the figure in Section C-l of this appendix* The 
formula for F^_j- 

F 2-5 = F AD^ d S + d ■» d S* r A* r S ) " F 2-7^ d A ) 

where an estimate of F^^td^) is made using (1-1). 

7. View Factor F 

Since node 2 represents a concave surface it radiates to 
itself. The view factor F^_ 2 is found from the equation 


F 2-9 + 


2-1 


+ F 


2-2 


1 , 


and is given by 

F 2-2 = 1 “ 2F AD (d S ,d S’ r S’ r S^ * 

8. View Factor F 3 _^ 

Node 3 represents all the material surface on the front 
of plate S. The calculation of F 3 _ 4 takes several steps. First we 
find F q . from the equations ' 

7 ~4 
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F - 1 > F 

* 9-4 ~ 1 * 9-13' 

F 9-13' = F 9-13 + ^ F SD^ d ’ r S’ r A ,S o^ * 

The expression for includes only the six nearest holes surround- 

ing node 13. More distant groups of holes consist of 6 \J~3 S q , 6 at 2S q , 
12 at 1.5 n/T S q , etc. The view factors from node 9 to node 13 and 
to each of the above groups are shown in Fig. 1-1 as functions of the 
plate spacing d. The contributions to the over-all view factor Fg_^i 
from holes more distant than S q are negligibly small. Next, we define 
node 9' as the disks corresponding to the fronts of all the holes in plate 
S, and note that 


F = F 

* 9 -4 9- 


9-4 ’ 


The view factor F^ g, is then found from the reciprocity condition 


4-9' 


Z_ TT 

A . ^ 9 '-4 


The area ratio A^,/A^ is defined by 


l-f> 


where f c and f^ are the opened area fractions in plates S and A. 


Finally, the desired view factor F^_ ^ is found from the equations 


A 4 f 1-f S 

F 3-4 = Aj F 4-3 = (-Tf^ 


-1 


4-3 


F 4-3 = 1 


4 - 9 ' 
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SPACING BETWEEN ELECTRODES, d 


Fig. I- 1 . Terms in the over-all view factor F 





The resulting expression is 


3-4 " Vi -f 


1 -1^7 (1 “ F DD (d ’ r S’ r A> ' 6F SD (d ’ r S’ r A’ S o> 

A 


9. View Factors i» F 3 _ 7 > 

Nodes 5' and 7' correspond to the cylindrical wall and 
front disk surface of all of the plate A apertures. From the equation 


F 3- 13 ' " 1 ~ F 3 -4 " F 3-5 1 +F 


3-4 " 3-5 1 3-7' ’ 


we learn the sum (F 3 _ 5 , + F 3 _ ?I ). No good approximation has been 
found to date for either; the difficulty is apparently that the shielding 
due to plate A enters significantly. Until more time is spent on their 
calculation, we will take 


F 3-5' “ F 3-7' - 2 F 3- . 


J 0 - F 3-4> • 


View Factors Fj-_^, F,._g 

The view factors F $ F,-_ 5 are given by the formulas 


F 5-7 = F AD (d A' d A- r A’ r A ) 


F 5-5 " 1 2F 5-7 * 
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APPENDIX II - A MODEL FOR THERMAL BUCKLING OF THE 
ELECTRODES 

The simplified model of thermal buckling considered here was 
developed for the purpose of explaining the computer results for elec- 
trode buckling. An analytical solution of the model problem is obtained 
and shown to agree qualitatively with! the computer results. In the mod- 
el problem a peak also occurs in the curve of axial displacement versus 
initial axial position. 

A. Computer Results for the Electrode 

The buckling of .an initially dished 30 cm molybdenum electrode 
(0. 030 in. thick) due to a radial temperature gradient T(r) was calcu- 
lated (by a "finite element method" computer program) for various 
values of initial eccentricity e (see Fig. II- 1). The electrode was as- 
sumed to be simply supported at the edge and to have no temperature 
variation through the thickness. The computer program actually ana- 
lyzed a solid plate with one small (0. 075 in. ) central hole, but employed 
modified elastic constants E and v to reflect that the electrode is 
perforated. A graph of the calculated deflection q/d versus e/d is 
shown in Fig. II- 2. It is the negative slope portion of this curve for 
which an explanation is sought. 

B. Model for Thermal Buckling 

The model for thermal buckling consists of a triangle of rods, 
hinged at the vertices a, b, and c (Fig. II-3). Each rod has cross 
section A, Young's modulus E, and thermal expansion coefficient a. 

In the un loaded state (Fig. II-3(a)), the rods are in mechanical equilib- 
rium with no internal stresses. The loaded state (Fig. II-3(b)), of the 
model is reached by heating rods ac and be to temperature T, and 
in turn compressing spring S. Rod ab remains at zero temperature. 






The analogy of the state of Fig. II-3(b) with the deformed state of 
the electrode is as follows: (see Fig. II- 1) 

1. Rod ab provides the radial clamping of the outer, rela- 
tively cold region 1 on the electrode. 

2. Rods ac and be correspond to the thermally expanding 
central portion 2 of the electrode. 

3. The spring provides the bending resistance of the cap 3. 

The solution of the thermal buckling problem under the assumptions 
(e/d) 2 « 1, (^d) 2 « 1, (r\/d) Z « 1, aT « 1, p « 1, a T < 2 |3 


Vd = (e / d)aT 


(e/df + 2(3 - aT 


(II- 1) 


P = 


k/2 

EA/d 


The graph of this solution is shown in Fig. II-4 for 


a =10 6 /°F, 
T = 100°F, 

7 

E = 10 psi, 

A = 0.1 in. ^ 
d = 10 in. , 

k = 100 Ib/in. 


Upon examination of the force balance equations for equilibrium 
a physical explanation of the peak is discovered. Because of the spring 
the displacement t| is proportional to the axial component of the com- 
pressive force F c in rods ac and be. The axial component is 
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F sin ^ abc and sin £ abc is proportional to e. It turns out that F 

C Z - 1 2 C 

is proportional to (e + constant) and therefore r\ a e/(e + constant). In 

the next paragraph the exact fourth degree algebraic equation for rj is de- 

rived and solved for (r|/d) « 1. Additional assumptions of (e/d) « 1, 

aT « 1, |3 « 1 are made to further simplify the solution and obtain (II- 1). 

These assumptions apply to realistic situations as in the example of Fig. II-4. 


C. Derivation of Equation (II-l) 

The length L q of rods ac and be before loading are, from 
Fig. II-3(a), 


L, 

o 


fd + 


2 . 1/2 

e } 


(II-2) 


After loading by thermal expansion and the spring force the length is, 
on the basis of Fig. II-3(b), 


Lr 


{(d + i) 2 


+ (e + -n) 2 



(H-3) 


The compressive force F c in rods ac and be is equal to the force 
required to compress each from length L q (1 + a T) to length L and is 


F 

c 


L q ( 1 + aT ) - L, 

L (1 + aT) 
o 


(II - 4) 


The tension force F in rod ab is equal to the force required to stretch 
it from length 2d to length 2(d + |) and is 


F 


t 


EA 


2(d + £) - 2d 
2d 



(II-5) 


Two equations for the unknown displacements £ and r| are provided by 
equating the vertical forces acting at hinge a and the horizontal forces 
at hinge c. These equations are 
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d + £ F • = F 
L. c t 


(II-6) 


^ F C = 


(ir-7) 


A single equation for t| can be derived from (II- 6) and (-II-7) as follows. 
After substitution of the expression (II-4) for F^ eq. (II-7) is solved 
for L to give 


L = L (1 + aT) 

o e + A 




(II-8) 


where 


A = 1 + 




( 1 + aT). 


The relationship 


t _ 3 Ti d 3 _ _kZ|_ 

“ P e + t, (1-P) ’ P ' EATd » 


(II- 9) 


is found by dividing (II- 6) by (11^7), is then used to eliminate £ in (II-3) 
and obtain 


L = e +Vl-B) {d 2 + [e +T,d - MJ 2 } / 

Finally, equating (II- 8) and (11-10) yields 

L o (1 + aT) 6 e ^qA 3 " = ^ + [ 6 + 71(1-13)1 2 1 


(II- 10) 


1/2 


(II- 11) 


If the variables e and -q are normalized and r) ' is assumed 
small, i. e. , 


e' = (e/L. Q ), t }' = (rj/L o ) « 1, 

eq. ( II— 11) simplifies to 
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. — ' 


6 


2 1/2 

(1+ aT) {e 1 +y (l-(3) } = {1 +ti' {2e(l-p) + n' (1-P) }} (e'+ n 'A) 

= {1 +V {e 1 (1-p) (1-p) 2 }} (e 1 + y A) 

= e' + T1 > {A + e' 2 (1-p) } + . . . (11-12) 


Solving (11-12) for r|' gives 


„ 1 - 


e' aT 


e ,2 (l-p) + p(l + ip) (1 + aT) - a T 


Under the additional assumptions that 


e ' 2 « 1, aT « 1, p « 1, 


eq. (11-13) reduces to 


t\ 


t _ 


e 1 aT 

e' 2 + 2p - aT 


(11-13) 
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